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The field of the invention relates generally to systems and methods for estimating credit risk 
associated with credit instruments and portfolios of credit instruments. 

BACKGROUND OF THE INVENTION 

Credit risk models are used to estimate the risk associated with credit instruments and 
portfolios of credit instruments, such as bonds and loans. Various models have evolved 
overtime. Many run simulations to generate possible valuation scenarios. Simulations are 
typically generated using a Monte Cario method or variant. 

Credit risk models typically generate scenarios using probability distributions based on 
historical data. These probability distributions may be obtained in the form of a transition 
matrix from various sources, for example Standard & Poors™ and Moody's™ Such 
probability distributions are often referred to as (P), a real measure probability distribution. 

In credit risk models, transition matrices under (P) probability measure, or (P) transition 
matrices, are used to predict transition, or migration, from one credit class to another, "credit 
migration", over a given time period. A set of probabilities of migration, or transition, for 
different credit classes forms a transition matrix. 

Models are typically implemented in computer software as part of a credit risk system usmg 
compatible hardware. 
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Alternative matrices, models, syst( ms and methods for estimatmg credit risk and credit 
migration associated with credit in jtruments are desirable. 

SUMMARY OF THE INVENTION 

In a first aspect the invention provides a method of creating an adjusted (P) transition matrix 
of probabilities of credit migration for different credit classes. The method uses the steps of 

• obtaining a transition matrix under probability measure (P) 

• adjusting the (P) transition matrix to be consistent with default probabilities under 
martingale measure (0 computed, for instance, from bond market prices by making a 
column of default probabilities in the (/') transition matrix consistent with default 
probabilities calculated under (0 

• scaling the other entries in the (P) transition matrix to compensate for the adjustment 
while retaining the relative weights among non-default classes from the (P) transition 
matrix in the adjusted transition matrix. 

The method may adjust the (P) transition matrix in accordance with the following: 
= {/^i Pk } - are all credit classes 



- corresponds to default state 



= P../',.', J 



- adjusted (P) transition matrix for the time interval ] 



P./',.', J 



transition probability from credit class /< to credit class R 




P>:Jio.',) = 0 if J<K and Z^^-.^- (/(,,/,)= I 



SUBSTITUTE SHEET (RULE 26) 



3 PCT/CAOO/01440 

^0=Q(C) = ^g(M^,..Z,) 

f^"^ - obtained (P) probability of default for the credit class R 

(from the obtained (P) transition matrix) 

- (0 probability of default for the credit class R 



G 



- a set of credit classes (can be the whole set R = R^] 

or some subset 



1^1 - denotes the number of elements in (} ■ if Q is empty, put 

Co=0 

^( ) - standard normal cumulative density function. 

Alternatively, the method may adjust the (/') transition matrix in accordance with the 
following: 

=\\P".,VvJj\\ - adjusted (/>) transition matrix for the time interval [z^./,] for 
obligor n 

P'/'oA) - obligor's n transition probability to migrate from credit class 

R, to credit class R^ 

\n{z: + c;\ i = k 

Pk.j (toJ^) = 0 if J < K and pl ^ (t^,t^ ) = I 
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- obtained {P) probability of default for the obligor n over 
interval [/q,/,] if it is in the credit class at time /„ 

- ^C^) probability of default for the obligor over interval 
[/„,/,] if it is in the credit class R^ at time 

''' - instantaneous risk-free interest rate at time / 

y"< - obligor n stock log-return over interval [/q./J 

- obligor n drift function of time / 

In a second aspect the invention provides a credit risk model for the risk analysis of complex 
portfolios of credit instruments. The adjusted {P) transition matrix may be used to value 
credit swaps and other complex credit instmments The model may have a simulation engine 
to generate valuation scenarios based on the adjusted (P) transition matrix The simulation 
engine may generate scenarios based on a Monte Carlo method. The model may be 
implemented in computer software on a computer readable medium. 

In a third aspect the invention provides a method of pricing a credit swap. The method 
utilizes a number of steps: 
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. {P) transition matrix over a risk horizon [/„,/,] and the {Q) default probabilities over the 

same period [/„,/,] are obtained 
• an adjusted (P) transition matrix and (0 transition matrix over are created 

. a first series of forward (0 transition matrices for standard observation times 
t-i„t,.t,,... is built from the (0 transition matrix and an initial term structure of 
default probabilities under O over [/„,/,] 
. the first series .s translated to produce a second series, having premium payment dates 
as the observation times by working wth 'fractional' powers of transition matrices for 
subintervals 

. a third series of transition matrices is derived from the second series for co-migration on 

credit classes for pairs of obligors 
• a risk-free discount factor and probabilities contained in the third series are used to 

discount cash flow for each payment date over the period, and the discounted cash 

flows are summed for all payment dates. The sum is subtracted from a payofT function 

for fixed payments to give the credit swap price. 

In a fourth aspect the invention provides an extended credit risk model for the risk analysis of 
complex portfolios of credit instruments with stochastic interest and foreign exchange rates. 
The model may have a Monte Carlo simulation engine to generate valuation scenarios based 
on the obtained or adjusted (/>) transition matrix The simulation engine may generate 
scenarios from a joint lognormal distribution of the following risk factors: 

• country/industry market indices 

• obligors' idiosyncratic components 



SUBSTITUTE SHEET (RULE 26) 



6 PCT/CAOO/01440 

• base forward rates j ^^^^^^^^ 

corresponding to given maturities /„./ 

. cred.t spread .v,(0 for the h.ghest credit class and for matunty equal to the first 
bucket point i. on the corporate zero curve for each currency 

• incremental spreads A.y(/ )= v(/ )- V (i\ i-~> k' \ p 

^ ^.VJ ■\,l',h/-2 A -1, for maturity A (by credit class 

and currency) 

. .l,e difference i,„ ,(,).«,. ,(,)->, ,(-,) be,„ee„ ,he spreads for ,he Icwc, credi, 
rating for the longest and shortest maturity (by currency) 

• foreign exchange rates. 

From risk factors logreturn.s for each simulation run the model may compute: 

• obligors- credit ratings at risk horizon based on country/industry market index logreturns 
and idiosyncratic components 

• base zero interest rates for all currencies based on forward rates 
/(/„.0./(',./:)... ..,/(/, ,./J 

. credit spreads .v, (,, ) = .v ,(,,). ,Vv ) for maturity . , = 2....A' - 1 (by currency) 

. credit spread .v, = ^^.^^ ^^^^^^^^ 

maturity (by currency) 
• the rest of the spreads for maturity r. by partitioning the interval [o..s,,Oj] 

proponionaiiy to .v (/, ) by taking s/,^ j = ,j ,^ j!LeJ^ (by currency) 

. spreads for other maturities t, , by linear interpolation between .v.(/,) and .v(/J 

(by currency) 
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. corporate zero cun.es at risk h.rizon by aading the appropriate spread to the base curves 

(by currency) 
• foreign exchange rates. 

In other aspects the invention provides various models, matrices, software, systems and 
methods incorporatmg combinations and subsets of the aspects described above, or utilizing 
the aspects described above. 

BRIEF DESCRIPTION OF THE DRAWINGS 

For a better understanding of the present invention and to show more clearly how it may be 
carried into effect, reference will now be made, by way of example, to the accompanying 
drawings which show the preferred embodiment of the present invention and in which: 

Fig. I is a block diagram of the structure of a credit risk model according to 

the preferred embodiment of the invention. 
Fig. 2 is an example of Portfolio Summary (analytic engine was not used) 

output from the credit risk model of Figure I, 
Fig. 3 is an example of Obligor Summary output from the credit risk model 
of Figure I, 

Fig. 4 is an example of Distribution of the Portfolio Forward Value output 

from the credit risk model of Figure I, and 
Fig. 5 is an example of Marginal Risk vs. Exposure output from the credit 

risk model of Figure I, and 
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Fig. 6 is a block diagram of the structure of a credit risk model with 
stochastic interest and foreign exchange rates according to the 
preferred embodiment of the invention 



DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Probability distribution (probability measure) (P) is a probability distribution of risk factors 
based on historical data. In credit risk models, probability distribution (P) is used to predict 
transition, or migration, from one credit class to another, ''credit migration" A group of 
probabilities of migration, or transition, for different credit classes forms a transition matrix 
A series of transition matrices can be built for different time intervals and for different 
obligors. 

In the preferred embodiment, a given, or obtained, transition matrix under the probability 
measure (P) is adjusted to be consistent with default probabilities under maningale 
probability measure (O). In the preferred embodiment, the (P) transition matrix is adjusted 
by making a column of default probabilities consistent with default probabilities calculated 
under (O) The other entries in the transu.on matr.x are scaled to compensate for the 
adjustment while retaining the relative weights among the non-default classes in the transition 



matrix. 



The adjusted (P) transition matrix may be incorporated into a credit risk model to analyze 
complex portfolios of credit instruments, including credit derivatives. For example, the 
model may use the adjusted (P) transition matrix in place of the given (P) transition matrix in 
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a simulation engine to generate scenarios. The simulation engine may generate scenarios 
based on any Monte Carlo method. 

The model may be implemented in computer software as part of a credit risk system usmg 
compatible hardware. 

The adjusted (/') transition matrix may be used in the pncing of credit swaps. Credit swap 
pricing may be incorporated into the credit risk model as it is consistent with the 
methodology used in the model. 

This description has four Parts. 

Part 1 describes assignee s basic credit risk model. Credit VaR'^^', which utilizes deterministic 
interest and foreign exchange rates. Each obligor s standardized equity returns are 
decomposed mto weighted average of market indices returns (multi-beta model) where the 
weights are specified to appropriately reflect the obligor's participation in the corresponding 
markets and to model obligor s idiosyncratic returns. Part 1 has Appendix 1. 2 and 3. 

Part 2 descnbes the derivation of an adjusted (/>) transition matrix and related matrices. In 
particular the relationship between P and O probabilities for obligors to migrate from one 
credit class to another including probabilities of default is derived. An adjusted transition 
matrix and related matrices are created based on the derived relationship The adjusted (/>) 
transition matrix and related matrices allow portfolios of different instruments, including 
credit derivatives, to be analyzed. In the preferred .mboc.iment, the adjusted (P) transition 
matrix and related matrices have been incorporated into an extended version of the credit risk 
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model of Pan 1. Part 2 also describes a method of pricing credit swaps using the adjusted {P) 
transition matnx and related matrices. In the preferred embodiment, the method has been 
incorporated mto the credit risk model of Pan I The adjusted transition {P) matnx, related 
matrices and correlations between equity returns are used to compute joint probability 
distribution of obligors" credit migrations. 

Pan 3 describes an extended version of the credit risk model that allows for stochastic interest 
and foreign exchange rates. The model generates correlated interest rates, foreign exchange 
rates and obligors' credit migrations according to their joint probab.l.tv distribution. As 
stochastic rates are permitted valuation must be performed for each mstrumcnt m a generated 
scenario 



It will be recognized that the adjusted {P) transition matrix and related matrices, pricing 
method and stochastic features may also be incorporated into other compatible credit risk 
models. 
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Part 1 - Basic CreditVaR Methadolosy 
1. Introduction 

CreditVaR'^' is assignee s proprietary model for measuring and analyzing credit risk in a 
portfolio context. This description is made in reference to CreditVaR as the preferred 
embodiment of the credit risk model. Those skilled in the art will recognise that many of the 
principles described herein are equallv applicable to other credit risk models, such as J.P. 
Morgan's implementation of CreditMetrics™. 

The methodology used in CreditVaR is based on credit miuraiion analysis and on the 
assumption that an obligor's credit migration is driven by its asset value For additional 
discussion of credit migration analysis, please see CreditMetrics''^' - Technical Document, 
April 2, 1997, J P Morgan & Co Incorporated that is hereby incorporated by reference and 
available at www creditmctrics.com 

The meihodology siruciure and implementauon of CreditVaR arc shown in Figure 1 
Information required for credit risk analysis includes: 

Obligor data: 

Information about obligors is organized into a database containing details of their 
credit ratings, industries, and countries. 

Portfolio data: 
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Information about financial positions is organized into portfolios of exposures. It 
allows to cover different types of instruments such as fixed income instruments, loans, 
commitments, letters of credit, etc. 

Market data, transition prohahilitics and corrclatwiis: 

These include yield curves, spread curves, foreign exchange rates, transition 
probabilities from one credit rating to another for different credit rating systems 
(Moody 8 states. Moody 18 states. S&P 8 states), correlations bctw/een market 
indices 

The CreditVaR methodology may be implemented both analytically and as a Monte Carlo 
simulation, [t calculates two risk measures: standard deviation (in the analvtic and simulation 
engines) and percentile level (in the simulation engine) of the portfolio value distribution for 
a given time horizon The regulatory capital calculation is based on the percentile level, 
therefore, only simulation engine is used for the capital calculation Accordinuly. in the 
preferred embodiment the analytic engine was removed. Following sections describe the 
methodology and a possible implementation in more details 

2. Risk Statistics and iVIesisureinent 

The CreditVaR. methodology can assess the impact of changes in debt value for a given 
portfolio due to credit quality movements of issuers (obligors) - including downgrades, 
upgrades and possible defaults - which can occur within the time horizon, typically a period 
of one year. 
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There are two computational engines developed in CreditVaR; Analytic Engine (AE) and 
Simulation Engine (SE), Both are capable of producing mean and standard Ucviaiion for the 
portfolio s forward value (risk horizon later). The two sets of numbers should be quite 
comparable as the same modelling framework is embodied in both engines. Though 
distribution of value for a ponfolio of bonds is typically very much skewed, these statistics 
contribute to some understanding of credit risk present in the ponfolio. 

The Simulation Engine generates a distribution of forward ponfolio values, and therefore has 
the advantage of estimating anv statistic, including ones deemed useful in risk measurement 
In panicular. Vulue-ai-Ri.sk for any percentile level, together with a confidence band, can be 
computed. (Another risk statistic that can be produced from SE is average shortfall which 
is the expected loss given that losses exceed a given level. The expected excession of a 
percentile level, which is the expected loss given that the loss is more extreme than the given 
percentile level, is another very useful risk measurement that can be obtained using SE 
engine ) Results become more accurate as number of simulations increases, but will exhibit 
cost in terms of time and memory as with any Monte Carlo method. 

With additional computation, both engines AE and SE produce mar}rinal standard deviation 
by each obligor, i.e., the difference between the standard deviation of the entire ponfolio and 
that of the ponfolio ex the obligor. Plotting marginal standard deviation expressed as 
percentage of mean value for the given asset against the exposure size for each obligor, it can 
be used to identify that pan of the ponfolio which has concentration of risk, large exposure 
size and high percentage of marginal standard deviation. All of these measurements are 
usefijl in risk monitoring, as well as for capital allocaiion purpose. 
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The current CreditVaR method can be extended to produce measurement of risk which are 
due to obligor-specific returns. Risk due to market indices, in fact to any subset of market 
indices, can also be obtained. 



3. Pre-processing 



The methodology requires that the three types of data mentioned m Introduction be given as 
inputs. They are then used to prepare the portfolio in an input format applicable to the 
CreditVaR Engine(s). We now indicate the required pre-processing 

For each instrument in the portfolio, there should be sufficient information to calculate its 
forward value for any credit migration with respect to a credit rating system chosen for 
CreditVaR Specifications such as coupon rate, term to maturity, yield curve, spreads, etc. 
are needed Recovery rate and its standard deviation, i e, the seniority of the debt, are also 
required in case default occurs. We assume that obligor's exposures can be aggregated or 
netted across instruments having same seniority 

The credit quality movements can be tabulated in a transition matrix (for the given risk 
horizon) of probabilities for an obligor in one rating to end up in other ratings, including 
default, or remain unchanged in its rating. An additional ingredient essential in the pre- 
processing is a module which quantitatively measures the ponfolio effect of credit by 
accounting for correlations of asset returns for the obligors. Once equipped with these asset 
correlations (which are on average typically between 20% and 35%), joint rating changes can 
be modelled across these obligors. 
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Ultimately, it turns out that in t lis model the only parameters which afTect the risk of a 
portfolio are the two already men.ioned. namely, the transition probabilities for each obligor 
and the correlations between asse; returns. In this framework, the formulation of risk 
estimation is then reduced to considering the '"standardized- (mean 0, standard deviation 1) 
asset returns, where the only parameters to be estimated are these correlations. 

It should be noted that interest rates and spreads are assumed to be deterministic in this Part 
1 In Parts 2 and 3 interest rates and foreign exchange rates are allowed to be stochastic As 
stochastic rates are permitted valuation must be performed for each instmment in a generated 
scenario. Part 3 describes a method and means for valuing such scenarios. 

As asset values arc in general difficult to observe, we use instead the correlation between 
equity returns as its proxy The obvious advantage of this is that equity data are more readily 
available than data of credit spreads or of actual joint rating changes. In order to keep the 
size of correlation matrix reasonably small to be useful in computation as well, the 
methodology relies on some form of component analysis. Basically, each obligor s 
standardised asset returns are decomposed into a weighted average of some benchmark" 
country-industry indices where the weights are specified to appropriately reflect the obligor's 
participation in the corresponding markets, plus a residual component that is associated to 
obligor-s idiosyncratic returns. These weights arc specified for each obligor as pan of input 
stored in Obligor database, (An example of such specification would be: Company X 
participates 70% m German Chemicals and 30% m German Electronics, and 15% of the 
movements in X's equity returns are company-specific) In this way, the methodology 
incorporates the mix in obligor s market capitalization as well as obligor-specific movements 
in returns when computing credit risk of the portfolio. Given these weights and a correlation 
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matrix for the indices, one can easily obtain a matrix of correlations between all of the 
obligors present in the portfolio. Appendix 2 presents in detail this methodology. 

We remark that in applying CreditVaR to a number of portfolios there arc several 
computational modules which can be made external to the pre-processing of individual 
portfolios. For example, the asset correlations between obligors can be stored in a database, 
and each time a portfolio is encountered for credit risk computation selectmg appropriate data 
can be considered as part of the pre-processing. 

4. Analytic tiigine (AE) 

Mean for the portfolio's forward value can be calculated in a straightforward way, using 
transition probabilities and forward bond prices which have already been valued in the pre- 
process Computing standard deviation requires more work, but it reduces to computing 
standard deviations of various obligor-pair subportfolios. The latter is straightforward as well 
since joint probabilities of credit quality co-movements arc made available in the pre-process. 
Marginal standard deviations can be computed in a similar fashion. 

5. Simulation Engine (SE) 

First, threshold levels of standardised asset returns representing credit rating changes are 
determined for each obligor in the portfolio using the transition probabilities, this part 
depends only on the obligor's credit class and not on the obligor itself Using the asset 
correlations as the covariance matrix, samples of standardised multivariate normal random 
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variables are generated. These are then bucketed using the threshold levels, thereby 
generating scenarios of credit rating co-movements for the obligors. 

For each of these scenarios, portfolio is revalued quickly using pre-processed prices of 
exposures Each time a default occurs in scenario, i.e., when the sampled asset return value 
of an obligor is below the default threshold level, a random recovery rate is generated 
according to a beta-distribution whose defining parameters are governed by seniority 
associated to the obligor s bonds Finally, we obtain a distribution of portfolio values and 
from it the relevant risk statistics. 

6. Implementation 

A CreditVaR software program utilizes the methodology described above. The core 
computational engine is implemented in C++™. The interface for Windows'"^' could be a 
combination of Microsoft Access^^'/Microsoft Excel™ Windows version will consist of a 
Microsoft Access/Excel application and a DLL with the compiled computational engine. The 
selection of hardware is at the discretion of the user provided it has sufilcient processing 
capability and is compatible with the programming environment chosen. Other 
programming environments could be used as would be evident to a person skilled in the art. 

CreditVaR program uses a country-industry covariance matrix to compute asset correlations 
between the obligors in the portfolio. Matrices may be contained on a computer readable 
medium, for example on a hard disk, floppy disk, optica! disk, or random access memory, as 
appropriate for the particular computing environment anc. stage of use. Matrices may be 
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contained in an appropriate signal for transmission, such as across the Internet or another 
computer network 

Computation engine consists of 2 parts - analytic and simulation. Analytic solution produces 
total mean and total and marginal standard deviation of every obligor. Simulation produces 
total mean and standard deviation of every obligor, as well as Value-at-Risk numbers for the 
whole ponfolio. For initial implementation time horizon is taken to be 1 year. The model can 
be expanded to handle other credit instruments like loans, loan commitments, etc. 

inputs 

The current implementation of the Global Analytics Credit VaR model evaluates the credit 
risk of a bond portfolio based on the following information: 
1 Obligor data 

• Current credit ratings of the issuers (obligors), company-specific risk of the obligors, 
capitalization data 

2. Market Data 

• Current zero-rate curves and current spread curves corresponding to different credit 
ratings and currencies 

• Forward ( 1 year from now) zero-rate curves and forward spread curves corresponding 
to different credit ratings and currencies 

• Transition probabilities matrix for credit rating migrations 

• Covariances between the market indices underlying obligors' capitalization 

3. Spot and forward (1 year) exchange rates 

4. Market Data exposure data (PortEfolio) 
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• For bonds: principal amou it of the bond, coupon rate or spread, seniority (defined by 
recovery rate and its standard deviation), maturity of the bond, currency, etc. 

• For loans and loan commitments: total authorized amount, loan amount, letters of 
credit amount, coupon or spread, fees on letters of credit and unused portions, 
seniority, etc. 



Outputs 

The program produces the following outputs: 

1 Aggregated by obiigor: 

• Total present value of the bonds by a given obligor 

• Total expected I year forward value of the bonds (obtained analytically and by 

simulation) 

• Total standard deviation of I year forward value of the bonds (obtained analytically 
and by simulation) 

• Stand-alone VaR at different confidence levels 

• Marginal standard deviation and marginal risk of a given obligor 

• Delta standard deviation and delta VaR 

2. Totals for the portfolio: 

• Total present value of the portfolio 

• Total expected portfolio value in I year from now (obtained analytically and by 
simulation) 

• Total standard deviation of the portfolio value in I year from now 

• Skewness and kurtosis of the portfolio distribution 
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• Percentiles of the portfolio value (VaR numbers of the portfolio) and expected 
excession at 0,1%, 0.5%, 1%. 5% and 10% levels 

• 95% confidence intervals for the VaR values 

Marginal standard deviation and marginal VaR for a given obligor / are defined as 



and S7'J\,,^i_^^^^^^^^ , 
io without positions 



where S77\„^,, laR,,,,^, - are portfolio standard deviation and VaR 
^"K..., „«,,,„, - arc standard deviation and VaR of the portfoli 



corresponding to this obligor. 

Obligor's Delta standard deviation (Delta VaR) is defined as the sensitivity of the portfolio's 
standard deviation (VaR number) with respect to the aggregated obligor s , position value 
P„hhK«r, ^ multiplied by /J,,^,,,,,^ . 



DellciVcili,, 



rVall , 



DeliaSTD,,,,^^^^ and DcltaVuli,,,,,^,,^ satisfy following equations [3] 

In other words. DeUaSTD,„,^„^_ UMiaVaR,,,,^^ ) is an obligor , contribution into the 
portfolio's total standard deviation (VaR) 
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Appendix 1 to Part 1 

Sample Input and Output Data 

Obligor data table (example) - contains the following information: 

• Short and Full name of the obligor 

• Current credit rating of the obligor 

• Company specific % - the percent of the company's risk that is not explained by 
company's sensitivity to a particular country or sector of industry 

Country 1 (2, 3) and Industry I (2, 3) for each country and appropriate percentages - 
capitalization of the company 
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Name 




Rating 


Company- 
Specific % 


Countr 
y 1 


Capitalization 
in Country 1, % 


Industry 
1 


Capitalization 
in Industry 1, 
% 


AAAAKOT 


General Electric 


Aaa 


20 


US 


100 


GNRL 


100 


AAAABPD 


Hydro-Quebec 


Aa 


61 


CA 


100 


GNRL 


100 


AAAAAUD 


Ontano, 
Province Of 


Aa 


72 


CA 


100 


GNRL 


100 


28664 


Dharmala Saktl 
Sejahtera, PT 


Ba 


51 


ID 


100 


GNRL 


50 


29715 


Citicorp 


A 


53 


US 


100 


GNRL 


20 


29841 


Anwar Sierad, 
PT 


B 


51 


ID 


100 


GNRL 


100 


27508 


BOC Group Pic, 
The 


A 


52 


GB 


70 


GNRL 


50 


29842 


Dharmala Inti 
Utama 


Ba 


51 


ID 


100 j 


GNRL 


100 



Short 
Name 


Long 
Name 


Capital- 
ization 
in 

Industry 
1,% 


Indust 

0' 
2 


Capital- 
ization in 
Industry 
2,% 


Industry 
3 


Capital- 
ization 
in 

Industry 
3,% 


1 Count 

n/ 

2 


Capital- 
ization 
in 

Country 

2,% 


AAAAKOT 


General 
Electric 


100 

1 














AAAABPD 


Hydro- 
Quebec 


100 














AAAAAU 

D 


Ontano, 
Province 
Of 


100 














28664 


Dharmala 
Saktl 

Sejahfera, 
PT 


50 


AUTO 


20 


FOOD 


30 






29715 


Citicorp 


20 


BFIN 


80 










29841 


Anwar 
Sierad, PT 


100 


ELCS 


50 






US 




27508 


BOC 

Group Pic, 
The 


50 














29842 


Dharmala 
Infi Utama 


100 
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Current zero-rate curves for comp iting present values of the bonds f fragment ): 



Currency Maturity 
AUD 



— i ^ ' \^aa I Dd [ [-j | Qgg 



- ' v^. vy-^-^jvju I u.ujoooo I u.UDi Tbb U 067266 I 0 n777RR 



0.0578 63 I 0.056^ 0.057031 I 0.0588 45 I 0.065964 o T^fJ ^^^iriS^ 



0 06 2668 I 0.059883 I 0.060592 I 0 OB?0787o-07TWoQ7 i^frgT^ 
0.065801 I 0.062786 I 0.063617 I 0.0 657?n I 0 0 7.990 I TTTTiT^j—-^^^ 

n nRR/1£;9 n necoec r r\ r^rr-^r^r^ , ^ - ■- ^ ^ ouo 



0.068452 I 0.0 65865 I 0.066699 I 0.06959 7 I 0 079R4R I n nQnmi 
0.070039 10.068845 10"'-"" ' 



, ^-^^^w I 0 072824 I 0.0843 33 I 0 094867 

0 070246 I 0.07 7058 I 0.078239 I 0 081695 I 0 .095994 I 0 inR9«Q 

n nCOC 1 C i n n-lr-r'.n „ ] ! ^ ! "^"^ 



: ; -^.v^.^ ^^jj I u UOI030 I u.uao aaA I u 108289 

0.068616 1007 5739 I 0.077275 I 0.081163 I 0 .100634 I 0 1 1^7T7 

n nA A-rnn l n r\An-,r^r^ , „ . _ _ — _ ^ 



0.125788 



— — . ■ ■^■>/iv/.uuiioo|ui UUP J4 I U 1 1 1) / 1 2 

0 044700 I 0 0 40700 I 0.041400 I 0.04310 0 I 0.047600 I 0 n-^^ynn 
n n47Q77 I n nAioAc I n n. ',^ I ._: 1 — : ~ 



0.047977 I 0042845 I 0.043648 I 0.045248 I 0.051292 I 6"05i99 



— — : — . I v/.v/^o^Ho I u.uj i^a^ I 0.01)6991 

0 051217 I 0 0 45657 I 0 046252 I 0.048066 I 0.055 175 I n orT?^ 
0 055735 I 0.05 0484 I 0.051195 I 0 052688 I 0.061934 I n nM^iT 

n O'^^QAP.O I n ntiCCPT I n nee -HO . ,^ , ■ ~ 



0 070827 



0.059462 1005 5583 I 0.056418 I 0.058528 I 0.068Q8R I 0 0775^ 
n ncTrne i n nc^o^^ , >^ >^ >— ! ' '^~>'* 



0.063016 I 0 059309 [ 0 060148 I 0063046 i 0.07335 



— . , — I ^.uu^u.<o I u u/ jjD I u UB3559 0 1 3q'l7R 

0.066228 I 0.06-2^6 _IO06 3372 I 0.0663 95 I 0.077919 I 0 n T^ IJ^tTn^ 
° '-M65104Tm6620^^ 



; . - ■ . ^^^^^ I w .uu3.<.<o I u.UO^ / U UU94()ST n 1R7Q/1A 

^^625^0.069235 I 0 o7o76iTo 074^9BT777^^7T^^ ^ 
n n4finR4 l n m^iano I n nc^^o , ^ „^ ^ iQ'aHo 



-— — : . w.^. I u u,HD3D I u.ua^ 1^5 I U. I 09529 I 0 187946 

0.046084 I 0.051 908 I 0.052608 ! 0 054308 I 0.058808 I 0 06400^ i n n7..n. 
n DilflnRQ I n nt;i7Q7i n nn-^i^ n,^ , . _ ■- " •J'*uo 



-rrrrr—— -.^-.^^w^ , ^ ^^ ^^^u , u.uooouo I u ub49U8 0 075408 

' nun ! °°''''' JM ^P I 0.0541 9/ I 0.060242 I 0 0B592 1inn7°7.. 
U.051818lO0534091O05 4000To055812l0.0629T9T7n^fi 95^^ 
0.055428 I 0.055384 0.0560 91 I n 0.7-^7^ i n nAC7oc i n ...t/. " 'T-]^.^^ 



^^^^^^^ I , u.uooo 1^ I u.ub^91 9 0.069272 I 0 091916 

-^g^^gJ^g^^ilMS ^ 0.057574 I 0.06671570^7^9 54^^ 
0 0 5 8 7 7 1 070 5 7811 | 0.058638 I 0.060727 | 0,07018? 



— , , ^^^^^ I w.uuu I u.u/ u I 

0.062437 I 0,060102 I 0.060929 I 0 nfi37RR n 0710^^ 

n nR71CO n nt-At-r\r> i ^ ^ , 1 - I 



0.069212 



Z I v^ ww^/uu| U.Uf JJttJ 

0.064606 I 0,065561 I 0 068582 | 0.080094 



0.066522 I 0.06761 9 [ 0 070830 | 0.083975 



. w , I wuou U.UO03/3 

0.067655 I 0,069103 [ 0.072814 I 0.090877 



0,079251 I 0.120275 



0.083994 



0.090640 



0.095206 



— ; >^"i^w.^wi-Ti I I u. I U4o4o 

0.025700 [ 0,0 18100 1 0 018800 1 0.020500 1 0.025 000 1 nmTT?^ 
00 29251 I 0 021334 I 0,022137 I 0^^? ^73Rinn9Q77n I n n^g^^g 
0,029647 I 0,024503 I 0,0251 01 i 0.02691 3 I 0 033995 I 0 04Q35~ 



, v.,v.^v^wo I w,uz3 i ui I u.u^byii 0.0339 951 0.040351 

0,029943 I 0.030250 I 0.030960 i 0,032465 i 0.041641 I 0,04882T 



0.138775 



0.156078 



0.187307 



0.049285 



0.062777 



— . — , w,wo^->u^ I \j xm ic* I \J VWo^.! 0 079790 

^^j°iZ2!M3534 5 10^36^^ , ° ^gj^,^ 

0 03101? 0 04niBQ I n m-TT i n n--,^^^ ; „ - " 



— — — : . , , , w,uuo^o^ I Ob U 056809 I 0 0974R3 

0 031012 I 0,07 0T8nMl 0F^ 043866 I 0,0541 16 I 0 OR l ^^ft^TTg?^ 



. ■ -.- • ^ iw^i , w UHJOD D I U.UD41 lb U 064209 0 1 1R'iin 

[M31123 i. 0.04^i27T?^^ 9 , 0.048570 | 0.060025 I 0 - 0705,5,5 n -^..n" 
A231259 I 0 077376 I 0,048456 I 0,051bT 9 | 0.064593 I 0 075804 n .9? 
0 031582lT055377 I 0,05690oTonfimHI~^ ' ' 



0.043600 I 0 040200 I 0.040900 I 0042600 | 0.0471 00 



I 0.096276 1 0,16215 



' "."-■"-^'-'w I u.uH-iouu I ^J. v^ / 1 uu 

0 045953 0 041731 I 0.042534 I 0.044133 | 0.050175 



- . . , w. yyT^^o^ I u.uHH [ I u.ua ul /5 

0.050143 I 0.045279 I 0.045874 I 0.04768 9 I 0 0.54801 
0.054534 I 0 049986 I 0.050697 I 0.052191 ^'^ 



0 058249 I 0.054607 I 0.055441 I 0.057548 



0.061660 I 0 058466 I 0.059305 I 0.062197 



I .^^ , v/.v/^j^u^ I u.uo^ I a/ U .U/^494 

0.065009 I 0.062601 I 0.063563 I 0.066605 I O07R?n4 



0.066847 I 0.064869 I 0.065974 I 0.069215 I 0.082483" 



^. ^.y^yj^^ ^^ : u.uoja/'ti u.uoa-^i5| 0 0 82483 

0.070431 I 0.074380 I 0.076066 I 0.080398 I 0.102531 



0.083789 



0.099872 



0.117637 



0.138390 



0.154962 



.187930 
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( I year; zero-rate curves tor computing forward values of the bonds ftraat 



Currencv 


Maturitv 


Aac 


Aal Al Baal Ba 1 Bl Caa 


^ 


C 


0.05555S 


0.053202 1 0.054107 I 0 055601 1 0.063192 1 0 068470 1 0 085RS4 






0.05555S 


0.053202 1 0.054107 1 0.055601 1 0 063192 1 0 068470 1 0 085R'^4 







0.060177 


0.057528 1 0.058066 1 0.059937 | 0.068371 1 0 074863 1 0 103750 







0.065036 


0.061292 1 0.062003 1 0.063435 1 0.073895 1 0 081 370 1 0 1 18040 


— 


§ 


0.067908 


0.064212 1 0.065066 I 0.067236 1 0.077591 1 0 087239 1 0 1 33997 


— 


9 


0.070155 


0.067162 1 0 068010 1 0.071044 1 0.081944 1 0 092575 1 0 1 S3fiR3 


— ^yt 


li 


0.071248 


0 069894 1 0 070867 I 0.073988 1 0.086007 | 0.096861 1 0 1666R4 





19 


0.071 148 


0.078271 1 0.079478 1 0.0830291 0 097859 1 0 110491 1 0 165043 


!££ 


29 


0 069150 


0 076488 1 0 078053 1 0 082017 1 0.102020 1 0 117422 1 0 163480 





0 


0.051264 


0,044994 1 0.045901 1 0 047400 1 0.054998 1 0.060293 1 0 07749r" 


BEF 


L 


0 051264 


0 044994 1 0.045901 1 0 047400 1 0,054998 I 0 060293 1 0 077496 


BEF^ 


2 


0.054490 


0.048145 1 0.048686 1 0 050557 I 0,058983 | 0.065477 | 0 0942RR 




4 


0,05851 1 


0.052945 1 0.053657 I 0,055099 1 0,065547 I 0,073032 1 0 109R10 


BEF 


6 


0.061942 


0,058084 1 0.0589421 0,061 122 I 0,071542 1 0 081230 1 0 128113 


GEF 


9 


0.065071 


0 061397 1 0 062252 1 0 065285 1 0 076252 1 0 086928 1 0 148053 


BEF 


14 


0.067782 


I 0 063984 \ 0 0649591 0.0680791 0.080118 1 0 091007 1 0 160574 


BEF 


19 


0.069923 


10 066404 f 0 057530 1 0.070851 1 0.084581 i 0,096219 1 0 194843 


BEF 


29 


0 063169 


0,070233 1 0 0717931 0 0758031 0 0957651 0.111506 1 0 192460^ 


CAD 


0 


0.050057 


0 051687 1 0.052592 I 0.054087 1 0 061679 1 0 066959 1 0 084146 


CAD 


■1 


0 050057 


0.051687 1 0.052592 1 0.054087 | 0 061679 I 0.066959 1 0 08414B 


CAD 


2 


0.054696 


0.054160 1 0.054697 1 0.056566 1 0.064981 1 0.071461 1 0 100264 


CAD 


4 


0.057776 


0 056255 1 0 056963 1 0 058392 1 0 068789 | 0 076227 10 11 2638 


CAD 


6 


0.060901 


0.058799 1 0 059647 1 0.061801 i 0.072090 1 0 081660 1 0 127933" 


CAD 


9 


0 064270 


0.061017! 0 061858 1 0 064847 1 0,075538 1 0.0861361 0 14604? 


CAD 


14 


0.068690 


0,065519 1 0 066492 1 0 059509 1 0.081531 1 0.092501 1 0 152067 


2AD_ 


19 


0.070443 


0.067297 1 0,068415 1 0,071707 1 0.085315 1 0 095824 1 0 193509 


£^ 


29 


0,071351 


0 0682021 0.069677 1 0.073458 1 0.092000 i 0 105251 1 0 191357 


CHF 


0 


0.032815 


0.024579 1 0.025484 1 0 026985 | 0.034552 1 0 039871 | 0 057028 


9ijL 


L 


0.032815 


0 024579 1 0.025484 1 0.026985 1 0 034552 1 0 039871 1 0 057028 


£t!L 


2 


0.031626 


0.027719 1 0.028266 1 0.0301 35 | 0 038522 | 0 045007 1 0 073527 


2t!L 


4 


0.031006 


0.033310 1 0 034023 1 0.035479 1 0 045844 1 0 053305 1 0 089455 


CHF_ 


6 


0.031270 


0 0382471 0.0390991 0 0412531 0 051598 1 0 061 155 1 0 107084 


2tiL 


9 


0.031604 


0 042673 1 0.043527 1 0 046495 1 0 057401 i 0 067953 i 0 127399 


^tlL 


14 


0.03151 1 1 


0 046548 1 0.047518 1 0 050605 | 0 052572 1 0 073430 1 0 141238 




19 


0.031553 1 


0 048939 1 0 050041 1 0.053304 1 0.066718 1 0.078210 1 0 168899 


CHF 


29 


0.031785 1 


0 056686 1 0 058239 1 0 052325 1 0.082248 1 0.098595 i 0 166578 


DEM 


0 


0.048311 1 


0.043265 1 0 044171 1 0 045559 1 0.053260 1 0 058552 1 0 075739 


DEM 


1 


0 048311 


0,043265 1 0,044171 1 0,045559 | 0 053260 i 0.058552 1 0 075739 


DEM 


2 


0.053430 


0 047828 1 0.048370 1 0.050243 1 0.058573 1 0.0551 71 1 0 093976 


DEM 


4 


0.057285 


0.052447 1 0.053160 1 0 054503 1 0.055050 1 0 072535 1 0 109106 


DEM 


6 


0.060711 


0.057027 1 0-057885 1 0.060050 1 0.070451 1 0 080129 1 0 126889 


DEM 


9 


0.063686 


0 0605151 0.0613691 0.0643981 0.075353 1 0 086016 i 0 147006 


DEM 


14 


0.066555 


0.064220 1 0.065201 1 0.068341 | 0.080461 | 0 091432 1 0 161773 1 


DEM 


19 0.068085 


0.066183 1 0.067310 1 0.070635 1 0.084379 | 0 096031 1 0 194855 1 


DEMI 


291 


D.071368 


00755781 0.0772991 0.0817261 0.1044931 0 123304 1 0 192463 



forward cur\-es were obtained from the current curves by bootstrapping method. 
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Table of foreign exchange rates (spot and forward): 



Currency 


USD per Ccy, 


USD per 




spot 


Ccy fwd 1 r 


ATS 


0.0812843 


u.uo 1 1 y j4 


ADD 


0.7000000 


u.Do 1 ^\ 


BEF 


0.0277393 


u.u^ /oyi D 


CAD 


0.7058159 


0 6984693 


CHF 


0.7041758 


0 71 63023 


DEM 


0.5720824 


0 571 7753 


DKK 


0 1502765 


0 1495353 


ESP 


0.0067732 


0 0066965 


FIM 


0 1895950 


0 1 886958 


FRF 


0.1708701 


0. 1 706838 


GBP 


0.6898000 


1 .DJD OOUO 


HKD 


0.1293494 


u. 1 / yy 


IDR 


0.0002805 


U.UUUU3 / 1 


lEP 


1 .4936000 


1 A^c\^^fi~fr^ 
1 .*f 1 u / 0 / 0 


ITL 


0.0005838 


u.uuuo / 


JPY 


0 0079494 


0.0083655 


MXN 


0 1222046 


0.1 05581 2 


MYR 


0.2928258 


0.2089619 


NLG 


0.5076400 


0.5070257 


NOK 


0.1407856 


0.1371640 


NZD 


0.6266000 


0.5710057 


PTE 


0.0056032 


0.0055445 


SEK 


0.1317020 


0.1275255 


SGD 


0.6327912 


0.5657095 


THB 


0.0252048 


0.0166732 


USD 


1 .0000000 


1 .0000000 


XEU 


1.1342000 


1.1190554 


ZAR 


0.2069536 


0.1866202 
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matrix of country-industry indices (fragment, actual size 21 1 by 21 1 ): 



Covariance 
matrix 




Country 


AU 1 AU 


1 
1 

AU 


AU 


Index name 
MSCI Australia 
Index f.CIAU) 


Country 
AU 


Industry 
AU_GNRL 


AU GNRL 

0.000200 


1 AU BFIN 

3 


1 AU BMEC 


AU CSTR 


ASX Banks & 
Finance Index 
(.ABII) 


AU 


AU_BFirj 


0.000165£ 


0.000370 

S 






ASX Media 
Index (.AMEI) 


AU 


AU_BMED 


0.0002384 


0.000144 

4 


0.000561E 




ASX Building 
Materials Index 
(ABMI) 


AU 


AUCSTR 


0.000171 


8.437E05 


0.0001 92 


0.00042T6j 


ASX Chemicals 
Index r.ACII) 


AU 


AU_CHEM 


0 0002567 


0.000162 

4 


0.0003524 


0.00023 


ASX Energy 
Index ( AXEYl 


AU 


AU_ENRG 


0.0001518 


0 000107 

9 


0.0001871 


0,0001241 


ASX Food & 

Household 

Goods 

Index f.AFHI) 


AU 


AU_FOOD 


0.0002046 


0.000192 

3 


0.0002408 


0.0001694 


ASX Insurance 
Index ( AIII) 


AU 


AUJNSU 


0.0002008 


0.000195 


0.000256 


0.0001246 


ASX Paper & 
Packaging 
Index 
J.APPI) 


AU 


AU_PAPR 


0.0001811 


7.415E-05 


0.0002756 


0.0002049 


j ASX Transport 
Index (.ATin 


AU 


AU_TRAN 


0.0001754 


0.000138 

9 


0.00025 


0.0001775 


MSCI Austna 
Index (.CIAT) 


AT 


AT_GNRL 


5.911E-05 


9.668E-06 


0,0001002 


7.506E-05 


MSCI Belgium 
Index (.CIBE) 


BE 


BE_GNRL 


8.795E-05 


5.433E-05 


0.000145 


9.354E-05 


MSCI Canada 
Index r.CICA) 


OA 


CA_GNRL 


0.0001152 


0,000125 


0,000127 


7 762E-05 


Toronto SE 
Automobiles & 
Parts Index 
(TIPA) 


OA 


CA_AUTO 


7.868E-05 


8.753E-05 1 


8,144E-05 1 


4.66BE-05 


Toronto SE 
Financial 
Services 
Index r.TFS) 


OA 


CA_BFIN 


0.0001036 j 


0.000180 1 
5 


0.0001739 1 


0.0001015 


Toronto SE 
Index r.TCMB) 


OA 


CA_BMED 


4,429E-05 


1,239E-05 


4,918E-05 


6.247E-05 


Toronto SE 
Cement & 
Concrete 
Index ( TIPC) 


OA 


CA_CSTR 


0,0001121 


7.239E-05 


0.0001713 


8.249E-05 


Toronto SE 
Chemicals 
Index 

(.TIPZ) 1 


OA 


CA_CHEM 


0.000119 


J.000I42 

2 


0.0001386 


8.82E-05 1 
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Transition probabilities of credit u ting migrations (Moody-8 rating system): 



Rating 


Aaa 


Aa 


A 


Baa 


Ba 


B 


Caa 


Default 


Aaa 


0.9338 


00594 


0 0064 


0.0000 


0 0002 


0.0000 


0.0000 


0 0002 


Aa 


00161 


0.9055 


00746 


0.0026 


0.0009 


0.0001 


0 0000 


0 0002 


A 


0.0007 


00228 


0 9242 


00463 


0 0045 


0 0012 


0 0001 


0 0002 


Baa 


0.0005 


0 0026 


0.0551 


0 8848 


0 0476 


0 0071 


0.0008 


00015 


Ba 


0.0002 


00005 


00042 


0 0516 


0 8691 


0.0591 


0.0024 


00129 


B 


0.0000 


0 0004 


00013 


00054 


0.0635 


0 8422 


0019! 


0 0681 


Caa 


0 0000 


0 0000 


0 0000 


0.0062 


0 0205 


0 0408 


0 6919 


0 2406 


Default 


0 0000 


00000 


0 0000 


0.0000 


0.0000 


0.0000 


0.0000 


1.0000 



Note; Any transition probability matrix may be used (not necessarily 8-state). 
Example 1 

Sample bond portfolio 



Exposure Name 


Oblitjor Name 


Ratine 


Currency 


Recovery Rate 


912810EN4-11 


AAAAKOT 


Aaa 


USD 


40 


021654878440001 


AAAAAUD 


Aa 


CAD 


50 


19971027-00539 


AAAAKOT 


Aaa 


USD 


40 


027054867160001 


AAAABPD 


Aa 


CAD 


60 


027052668030001 


AAAABPD 


Aa 


CAD 


38 


021655365450002 


AAAAAUD 


Aa 


CAD 


32 


021655355640001 


AAAAAUD 


Aa 


CAD 


17 
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Exposure Name 


RecovRatc 
StdDcv 


Amount 


Maturity 


Rate 


Freq- 
uency 


912810EN4-11 


25 


5000000. 00 


15/11/2022 


7 . 625 


2 


021654878440001 


25 


15151000.00 


11/03/2003 


8 


2 


19971027-00539 


25 


14000000. 00 


15/08/2027 


6.375 


2 


027054867160001 




4 905000 . 00 


15/08/2020 


11 


2 


027052668030001 


23 


1300000.00 


15/10/2004 


11 


2 


021655365450002 


20 


25000000.00 


02/06/2027 


7.6 


2 


021655355640001 


11 


15420000.00 


01/12/2005 


8.25 


2 



"Rate" is the fixed coupon rate of the bond, "Frequency" is coupon payment frequency 

year Obligors 

AAAAKOT - General Electric 
AAAABPD- Hydro-Quebec 
AAAAAUD - Ontario, Province Of 

Sample output (obtained by running CreditVaR program on the above portfolio) 

Portfolio: Exp7 

Number of simulation runs: 100000 
Number of exposures: 7 
Number of obligors: 3 
Evaluation date: 14/04/98 
Started: 11:14:37 
Finished: 11:16:24 
Simulation time: 00:01:47 
Analytic solution time: 00:00:00 
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Portfolio Present Value = 71786803.00 



Portfolio Mean (analytic) = 70144493.80 
Portfolio StDev (analytic) = 523028.46 



Portfolio Mean (simulation) = 70143025.60 

Portfolio StDev (simulation) = 530826.76 

Portfolio Percentile (10%) = -105932.78 

Portfolio Percentile (5^) = -297306.87 

Portfolio Percentile ill) = -812778.79 

Portfolio Percentile {0.51) = -970228.53 

Portfolio Percentile {0.1%) = -4912661.85 



Obligor 


Rating 


Present Value 


Mean 
(analytic) 


Mean 

(simulation) 


AAAAKOT 


Aaa 


19190259.30 


19093143 . 19 


19093014 .25 


AAAAAUD 


Aa 


45880465.72 


44557820.87 


44556843.65 


AAAABPD 


Aa 


6716077.99 


6493529.74 


6493167.70 



Obligor 


StDev 
(analytic) 


StDev. 
(simulation) 


Marginal 
StDev 


Marginal 
Risk 


AAAAKOT 


73207.89 


73760.89 


3636. 56 


0 .019% 


AAAAAUD 


507005.01 


515008.37 


422137 .89 


0. 920% 


AAAABPD 


71014 . 55 


78555.92 


12079. 87 


0. 180-4 



Example 2 (actual program output) 
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Exposure data (fragment of input, actual size is 1863 exposures. 207 obligors) 



ID 


Exposure 
Name 


Obligor 
Short 
Name 


Credit 
Rating 


Currenc 


Recovery 
Rate 


Recovery 
StDev 


Amount 


Maturit 


i Rat£ 


Frequency 


15 


15656 


313122 


Aa 


AUD 


40 


25 


1 OOOOOC 


3/20/06 


6.4 




16 


15657 


313122 


Aa 


AUD 


40 


25 


1000000 


3/20/13 


6.57 
5 





17 


15658 


313122 


Aa 


AUD 


40 


25 


1000000 


3/20/1 8 


7 


2 


18 


15659 


313122 


Aa 


AUD 


40 


25 


1000000 


3/20/28 


7.05 


2 


19 


15660 


313123 


A 


AUD 


40 


25 


1000000 


3/20/99 


5,41 


2 


20 

IT 


15661 


313123 


A 


AUD 


40 


25 


1000000 


3/20/00 


5.41 




15662 


313123 


A 


AUD 


40 


25 


1 000000 


3/20/01 


5.64 




22 


15663 


313123 


A 


AUD 


40 


25 


1 000000 


3/20/03 


6 




23 


15664 


313123 


A 


AUD 


40 


25 


1000000 


3/20/05 


6.27 




24 


15665 


313123 


A 


AUD 


40 


25 


1 000000 


3/20/08 


6.53 


2 


25 


15666 


313123 


A 


AUD 


40 


25 


1000000 


3/20/13 


6.71 


2 


26 


15667 


313123 


A 


AUD 


40 


25 


1000000 


3/20/18 


7.14 




27 


15668 


313123 


A 


AUD 


40 


25 


1000000 


3/20/28 


7.2 




28 


15669 


313124 


Baa 


AUD 


40 


25 


1 000000 j 


3/20/99 


5.61 




29 


15670 


313124 


Baa 


AUD 


40 




000000 1 


3/20/00 


5.6 


2 



Example of the program output - Portfolio Summary (analytic engine was not used) 
(see Figure 2) 

Example of the program output - Obligor Summary 

(see Figure 3) 

Example of the program output - Distribution of the portfolio forward value 
(See Figure 4) 

Example of the program output - Marginal Risk vs. Exposure 
(See Figure 5) 
Appendix 2 to Part 1 
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Correlation model 
2,1 Market and specific components and country/industry weights 

Let us assume that company's standardized equity return r of the firm A can be expressed in 
the form: 

^ = /?u^u-^M, (1) 
where r^, is standardized return of the market component of r and r, is standardized return 
of the idiosyncratic component. and are independent. From this follows that 
coefficients p^, and satisfy condition 

fiu+fi^ = ^ (2) 
We also assume that return of the market component r„ is linear combination of 
standardized country/industry index returns r, , which are independent of r^. . 

r^,, =Pr,+... + P„r„. (3) 

cor{r,,r^)=0, i = ] „ 

Equations (I), (3) can be obtained by regressing A equity returns against index returns in a 
usual way. 



There is another way to represent r in terms of index returns r, and idiosyncratic return r. 
This representation is consistent with (l)-(3) but based on country/industry weighting 
coefficients. We show below how this representation can be obtained from (])-(3). 
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Denote by cr, cr, - standard deviations of the firm A equity return and index / return 
respectively, p,^ correlation between index returns r and (iJ-\ //) Then from (3) 

one has 



(4) 



where if, , w„ are given by 



* s 



(5) 



Coefficients 
percents. 



Vi'i + M', + . . . + = 1 , 

are called country/industry weights and can be expressed in 



From equation (4) follows that 



(6) 
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(7) 



: (8) 

Using (I), (4), (8) one can express standardized equity return in terms of country/industry 



/eights a 



Hence, to compute this representation one needs to have following parameters: 
country/industry weights, index standard deviations and correlations, specific coefficient fi. 
(coefficient for the market component can be computed from (2) as (i^, = ^\ - fi- ). in 
the current implementation CrcditVaR 1 derives parameters fi, and from "specific risk 
percenf .v defined as a fraction of total equity return changes explained by tirm-specific 
movements 



2.2 Correlation between two obligors 

Equations (I). (3) as well as equation (9) allow to compute correlations between equ: 
returns of any two firms. 

Let r", r' are standardized equity returns of two different firms B and C respectively and 
SUBSTITUTE SHEET (RULE 26) 



wo 01/52121 

34 
where 

^-"^A';^. A'=A';^.'=i....,// (12) 

Then correlation between r" and / ' is 

"4-"^'') = ZZA"A''p. (13) 

Example. 

Suppose we wish to compute correlation between two obligors B and C in the obligor table in 
Appendix 1 (Sample Input and Output Data). Obligors short names are "28664" and "29715" 
respectively. From the tabic we have 

Obligor B ("28664") 

Specific percent .v" = 0, 87 

Country/industry weights = 0.5 , v." = 0.2 , w'!^ ^^^^ = 0.3 

Obligor C ("297 15") 

Specific percent .v^ = 0.88 
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Country/industry weights w^j^ ^^^.^^ = 0.2 , 11^X8^- =08. 



Let us assume following standard deviations and correlations between market indices: 



Index 


Standard 
Deviation 


Correlations 


ID 
GNRL 


ID AUTO 


ID 
FOOD 


US GNRL 


US 
BFIN 


Indonesia, 
General 


003 


1 


04 


02 


0.01 


0.03 


Indonesia, 
Automobiles 


0.02 


0.4 


1 


0 1 


0.02 


004 


Indonesia, 
Food 


004 


0.2 


0.1 




001 


002 


US, General 


0,01 


0 01 


0 02 1 0 01 


1 1 0,3 


US, Banking 
and Finance 


0.02 


0.03 


004 


0.02 


0.3 1 1 



According to (7); 

CTj = 0.03' X 0.5' + 0.02' X 0.2' + 0.04' x 0.3' + 2 x 0.03 x 0.02 x 0.5 x 0.2 x 0.4 
+ 2x0.03x0.04x0.5x0.3x0.2 + 2x0.02x0.04x0.2x0.3x0.1 
= 0.000471 

o-g = a/0.000471 = 0.022 

CT^ = 0.01' X 0,2' + 0.02' xO.8' + 2 X 0.01 X 0.02 x 0.2x 0.8 x 0 3 
= 0.000279 

CTc = V0,000279 =0.017 

From(lO). (11), (12). 



J3^ = Vl-0 87' =0.49, 
= 0 87, 

.;^^^^ <o.v^.^.o.oM. ,o,, 0.5x0.03 ^^33 
o-g 0.022 

= fi^ <-^oroy..uTo ^ 0.49^1^ = 0.09, 
CTg 0.022 
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.0 49ililMl=0 27 
0.022 



P'-, = n/I-0.88' = 0.47 . 
= 0.88 , 



= K = 0.47^^ . 0.06 

CT^- 0.017 

= ^s.r.^rs... ^ 0.47MIM2 = 0.44 
o-c 0.017 

Finally, compute correlation between /•" and /•'' using (13): 

cor(r'\r'')= 0.33 x 0 06 x 0 0 1 + 0.33 x 0.44 x 0.03 + 0.09 x 0.06 x 0.02 + 0.09 x 0 44 x 0.04 
+ 2 X 0 27 X 0.06 X 0.0 1 + 2 X 0.27 x 0.44 x 0.02 
= 0.017 



2.3 Simulation of correlated equity returns 



Obligors correlation inatrix is used in the model to generate correlated samples of normally 
distributed random vectors of equity returns. In this section we describe the procedure based 
on Cholesky decomposition of the correlation matrix of index returns. 



Suppose there are A' obligors and their standardized equity returns are represented as (5 
(12) above): 

/^, =A'/', +--- + ^>„+/?.^'-/, k = \,...,K. (14) 

Denote by 
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Pzx Pzi 



p\ PJ 
p\ Pi 

Vp\ p^ 



P\n 
Pzn 



-i': II. ■ 



{ P\ 0 

0 Pi 



where (\, is correlation matrix between indices and I^. is A' x AT identity matrix. 



Then correlation matrix C of equity returns equals to: 

C = y-;(yTy7' ) = K{fV'Fr'fV)= W'E(rr' )V = ^V'C^iV (15) 
Let us assume that matrix H„ is obtained by Cholesky decomposition of the matrix C\, 

^-s,=H,,Hl,. (16) 

Denote by 

^ = ^'^.w (17) 
{H IS Kx{„ + K) matrix). Then 

C„ = H^Hl, 

C = W'H„Hl,W = HH' (18) 
are Cholesky decompositions of C\, and C respectively. 
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For simulation of a sample of K dimensional normal vectors /7 with zero vector of means 
and covariance matrix C (standard deviations of coordinates of R equal 1) we first simulate 
11 + K dimensional vector of standard independent normal variables 



■(a 



and then compute R as 

From (18) follows that R has distribution N{0,C): 
f-{R)=0, 

Equation (19) can be used for simulating correlated sample of equity returns in the model. 

Appendix 3 to Part 1 

Pricing Procedures 
3.1 Pricing Procedure for Floating Rate Notes 

Floating rate notes are designed so that the coupon payments depend on some current interest 
rate index (usually LIBOR or some other known rate), which is called reference rate This 
causes the value of a Floating Rate Note (FRN) to be close to par at all times. Credit risky 
FRNs pay coupons combined of reference rate plus seme spread over it, which depends on 
the credit rating of the issuer. 
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I' - price of note, 

/•, - expected future reference rate 

y, - spot long-term yields, used to discount the cash-flows at time / 
V - stated credit spread, fixed under terms of the FRN 
m - maturity of note in years 
// - coupon payment frequency 

Assuminy flat structure of credit spreads, the price of the note per 100 of face value, is 

Here, r, is forward reference rate for the period between /-th and (i IJ-st coupon payments. 
Assuming the reference rate is annually compounded, the way to compute forward reference 
rater, is the following: 



where R, is the rate on the reference cui-vc for time point /, . 

Within the current framework of CreditVaR I model, we will take existing current and 
forward corporate curves as a proxy for y„ since they incorporate market required spread. The 
reference rate R and credit spread .v should be supplied as parameters of the FRN. 
Note that for computing present value of an FRN we use the supplied reference curve, but to 
compute forward I -year value of an FRN, a forward reference curve has to be obtained. The 
rates r, in the computation of forward value of an FRN are therefore "forward forward" rates. 
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3.2 Pricing Procedure for Interest Rate Swaps 

Swaps can be characterized as the difference between two bonds, which can be either fixed or 
floating rate In the swap agreement, principal is not exchanged, which makes swap different 
from a bond from the point of view of credit risk Another difference is that in the event of 
default by the counterparty, the sign of the value of a swap determines weather the other 
party pays or receives some recovery amount (this depends on netting rules). 

Assumptions: 

1 ) There is no more than one credit event per year 

2) The credit event happens one year from evaluation date 

3) Netting is done for positions with the same obligor and seniority 

4) One recovery rate is applied to a netted position 

5) Implementation allows different recovery models depending on the sign of exposure to a 
bank, including recovery rate 0 if the value of the swap to the bank is negative 

6) A bank's credit events are not modelled fi e., the bank's credit rating is constant) 

The notation is as follows: 

V- value of the swap to the bank, 

Bi - value of the bond modelling the cash flow received by the bank, 

B] - value of the bond modelling the cash flow paid by the bank, 

r, - recovery rate applied if bank receives in the event of default by counterparty, 

r J - recovery rate applied if bank pays in the event of default by counterparty, 
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O - notional pincipal in swap agreement. 



The value of the swap to the bank then is: 

F = 8, - B, 

B, and Bz are present (or forward 1 year) values of the underlying bonds, obtained by 
standard procedure for fixed or floating rate bonds Ideally, swap bears credit risk only if the 
value of the swap to the financial institution is positive In practice, however, if counterparty 
defaults, it might not free the bank from liability to that counterparty (due to netting rules). 
Therefore, in case of default, the recovery amount is calculated as follows: 



If F 0, the recovered amount is ri*V, 
If K • 0, the amount paid is ri*K 



Here, we don't take into consideration the event of default or credit rating migration by the 
bank. Recovery rates r, and are determined from legal considerations at the time of default. 
We can use Beta distribution again to simulate recovery rates r, and (in current 
implementation, we take rr-r.). Alternatively, r, can be put to 0 if the swaps with negative 
exposure are excluded from the computation. 



3.3 Pricing Procedure for Loan Commitments 



A loan commitment is composed of a used portion and an unused portion. Used portion has a 

Loan (drawn flinds) and a Letter of Credit (undrawn ftjnds) parts. Interest (usually LIBOR 
plus spread) is paid on the Loan portion of a loan commitm :nt, LIBOR plus LC fee is paid on 
the Letter of Credit part (throughout its lifetime), and an Urused Fee is paid on the unused 
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ponion of the loan commitment. There is also a Facility Fee applied to the total loan 
commitment amount throughout its lifetime. A I -time Upfront Fee paid on the total loan 
commitment amount at the initiation can be ignored in the Credit VaR I framework since it 
bears no credit risk 

The amount drawn at the risk horizon is closely related to the credit rating of the obligor. If 
the obligor deteriorates, it is more likely to withdraw additional funds. On the other hand, if 
its prospects improve, it is unlikely to need the extra borrowings. 

Let 7 be the total (authorized) loan commitment amount; A - the amount drawn (in loans); LC 
- the amount in Letters of Credit, U- unused portion of the loan commitment. 
We have that 

T = L + IX: + U . 

Let / be LIBOR and .v be spread over LIBOR, so (I ■ s) is the interest paid on the loan; Ic - the 
fee paid on the letter of credit, ///- the fee paid on the unused portion;^- facility fee applied 
to the total authorized amount of loan commitment. Then present value (in the CreditVaR 
framework) is computed as follows we treat the loan and letter of credit parts as floating rate 
notes (because of the way interest is paid), facility fee and unused ponion fee are valued as 
bonds with no principal repayment. 



where is the annualized rate taken from the corresponding current corporate curve used to 

discount cash flows (depends on credit rating of the obligor in the CreditVaR framework); // 
is the frequency of coupon payments, and m is the number of years to maturity. Forward 
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value (FV) \s computed similarly, with .v being the fonvard (1 year from now) rates taken 
from the corresponding corporate ourve 

If downgrade or default happen to a particular obligor, it is likely to withdraw additional 
funds from the loan commitment line. To model this we assume that there is an "expected" 
drawdown (J-D), or "average commitment usage" if the obligor changes credit rating (see 
Credit Metrics Technical Document, p 45 for appropriate information) 



Credit rnting 


Averaue commitment usaue (ED) 


Aaa 


0.1% 


Aa 


1 ,6% 


A 


4.6% 


Baa 


20.0% 


Ba 


46.8% 


B 


63.7% 


Caa 


75.0% 



Then in the case of any credit event (downgrade, upgrade, default) the used portion of a loan 
commitment is adjusted in the following way: 



Downgrade 

Let X% {Y% ) is average commitment usage given in the Table above for the current (new) 
credit rating respectively, X <Y . Then new unused, loan and letter of credit portions are 

computed as 

l-X 

- new used amount is proportionally distributed between loan amount and 

letters of credit amount ZC^„ 
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Upgrade 

In this case X >V New unused. loan and letter of credit portions are computed 

' X 

= max (/.-((/„,^-(/),0) 



Default 

In the case of default it is assumed that used portion is the maximum of current used portion 
and 90% of the authorized amount 7" : 

^„.w = max(/, + /,C0.9x7-), 

Forward value of the loan commitment in the case of upgrade or downgrade is computed as 

In the case of default the loan commitment value is x where R is the recovery rate 
sampled from the Beta distribution (consistent with CreditVaR I framework), and total 
forward value is 
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Part 2 - CreditVaR Methodology: Credit Migration Process Under Two Probability 
Measures 



As described in Part 1, the CreditVaR is a model for measuring and analyzing credit risk in a 
portfolio context The methodology used is based on credit migration analysis and on the 
assumption that an obligor's credit migration is driven by its asset value. Correlations 
between equity returns are used to compute joint probability distribution of obligors' credit 
migrations. Each obligor's standardized equity returns are decomposed into weighted average 
of market indices returns (multi-beta model) where the weights are specified to appropriately 
reflect the obligor s participation in the corresponding markets and to model obligor's 
idiosyncratic returns. 

In the Monte Carlo simulation version the model calculates distribution of the portfolio 
values, statistical characteristics and risk measures defined as percentiles at given confidence 
levels (95%, 99%, etc). Simulation engine generates scenarios based on the "real world" 
probability distribution (P) of the risk factors. The distribution is assumed to be lognormal 
with parameters computed from the historical time series data. For each generated scenario 
the portfolio value is computed under the martingale probability measure (0 The 
relationship between these two probability measures is derived from the assumptions about 
stochastic processes for index and stock returns. In particular, we derive the relationship 
between P and O probabilities for obligors to migrate from one credit class to another 
including probabilities of default (transition matrices). This allows one to analyze portfolios 
of different instruments including credit derivatives. 
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In the CreditVaR framework it is assumed that credit migration process forms a discrete 
Markov chain with fixed moments of time f = {(/^ < {/. < ... < U, <...}, U., > 0 . In sections 
1-5 of this Part we present mathematical models of the migration process and define 
probabilities under P and (J measures to migrate from one credit class at time to another 
credit class at time ^Z^.,. Theoretical results of sections 1-5 are applied in sections 6-8 to 
compute P and O transition probabilities from market credit spreads and credit migration 
historical data. 

1. Stochastic processes for market index and stock returns under F measure 

Suppose that H = (h\B' H") is a standard /) -dimensional Brownian motion on a 

probability space {n,?,P) (B\H' B" are independent), where F = {/•;,/ >o} is a 

standard filtration of H . 

Let I = (/',/- 1"] be a set of all market indices Denote by V" the index value at 

time / > 0 . It is assumed that satisfies stochastic difTcrential equation (SDE) of the form: 

di T = r," (//,"' J/ + e"dB] + o;':dB- + . . . + e^,,dB',' ) , / > o 

(1.1) 

r„'" = v"', /n = l,2,...,A/, 
where 0" are piecewise continuous deterministic functions of /, / = 1,2 D . 



SUBSTITUTE SHEET (RULE 26) 



PCT/CAOO/01440 

Let us fix some time interval {/„/,] corresponding to two consecutive time steps from f : 
'o=U„ /, =U,.„!>0. Denote = ir(r,'")-in(F,:) be the index /" log-return for the time 
interval {/,,,/J, /„</</, Then from Ito's lemma and (1.1) follows that 

■ , (1.2) 
+ j0"dBl + d^^dB] +... + e'^adB° 

for /„ < / < /, with initial condition x,^ = 0 In a matrix form (1.2) can be written as 

-v, = a, + \d^dB^ ./„</</,, x, = 0 . (13) 

where 

xr^ 

dB^^[dB\ dB^) 

Let us consider obligors (issuers) S = (.V'.5= ^-). and let //; be an obligor stock 

price at time / We assume that //," satisfies SDE: 

dH", = H: (a^di + (7l4B] + al^dB] + . . . + o^^dB^ ) , / > o . 

(1.4) 

Hl^h\ /j=i,2 
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(1.5) 



where a", o-,", are piecewise continuous deterministic functions of time, / = 1,2,... ,Z) Then 
for the stock log-return = \n(H")- ln(//," ) one obtains 

+ f a",JBl + ct",.JB' +... + cT"DiJB'' , 

t^<t<t^, >'," = 0 . 
Similar to (1.3) this can be written as: 

y, =h, + jcT^JB^ , l,,<t<t^, y,^ = 0 , 



(1.6) 



where 

y, = {>'', >''T- 

h=(h: A;^)r 

2. Multi-beta model 

In the CreditVaR framework correlations between stock returns are derived from the index 
return correlations using multi-beta model. This is linear regression model that represents 
stock returns as linear combinations of index returns and stock specific (idiosyncratic) 
components. 

y^ = + /?;'x; + /3;-'x; + ... + /?,"■'" r," +< ,/„</</,, 

(2.1) 
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" = 1,2 A^, 

where < and .v; are uncorrelated: 6«v(C-v;')= 0 , and £(<)=0. = 1,2 A/, 

" = . In the matrix form; 

y, = A" + fi,x, + 6; , 

(2.2) 

where /?;'=(/?;•" P, =||A"H,.,, is a matrix of stock "betas", = 

a vector of idiosyncratic components which are uncorrelated with index returns: 
cov(f, , .r, ) = fcovt;. .r,"']!^ = 0 

(2.3) 

(here and below expected value E{) and covariance cov(.) are conditional to a-algebra 

It is assumed that rank{e,) = M<D for all (Note, that if w(.,)<,/,z? for 

some , = / then it can be shown that there exists time interval [r.rje [,,„,,]. / . [ij] and index 
subset I c I. 1.0, such that for all indices in I their log-returns can be represented as 
linear functions of index I.I log-returns w,th time dependent deterministic coefF.cients. 
Then index set I can be reduced to subset III on [,-,,=) ) Then from (1,8). (1.9) it follows 
that the {M x/W)-matrix efij is nonsingular on [/„./,] and 



(2.4) 

Indeed, from (1.3. (16) follows that 
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(2.5) 



^'"Kv,.x,)= \ avoids. 

Then (2.2) and (2.5) give ^^'^^ 
Therefore. 

]cT,e^d.s = p^\efild.s 

and (2.4) holds. 



components; 



(2.7) 



In the special case of constant volat.i.ty .atnces . and . one has 
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3. Credi. migration model under probability measure P 

as.„ed ,ba. a. ,ime ,;.f ea. obli.or .V-.S be,.,. .„ 

» ^"PPO- ".a. cedi, .i,a,io„ p.ces3 

forms a Markov chain on R with discretP t,m„c ^ . . 

d.screte t.mes T and absorbing state . Denote by R^^ , 

credit class of obligor .V at time {/ q;^ i 

t.me U e T S.m.lar to section 1, we fix a time mterval f/.,.,, j 

corresponding to two consecutive time steps from f , , , , 

/. 'i , / > 0 Suppose that 
K = f^ and ^,,) isaprobabilityor." to migrate to credit class R, at time 



(3.1) 



(3.2) 



Denote by 

a transition matrix of 5" for the time interval Then 

Zp:,{'..-^',)=i j = i K 

/^:-..('„->o=o, j = i fc-i. 

For each nondefault state /< and probabilities - 1 r H r . 

J-l,...,K define thresholds 

^ ^ ' ^ ' ' ^'^"'^^^d "o^al distribution as follows. 
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where /V(.) is standard normal cumulative density function. 



(3.3) 



- ^, IG; - /•-'(>',)) be a vector of standardized stock 
given by 



returns with correlation matrix 



where Z, , is a (yV x yV) diagonal matrix with stock 



(3.4) 



returns standard deviations 



VvarU') 0 

0 ... 



yvar(x^)^ 



(var(.) is a variance under measure conditional to o" -algebra ; 



(3.5) 



Similar to the CreditMetrics framework, the CreditVaR credit migration model is based on 
the assumption that credit class of the obligor r at time , is defined by the value of its 
normalized stock return n," and thresholds A," 
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Assumption: Given 7?-/^, - fc(/.,.O.AV,(/o,/,)J, where 



From the Assumption .t follows that joint migration process for obl.gors can be defined using 
a joint probability distribution of their normalized stock returns. Indeed, suppose that 

/e,;;=\, .VeS, ,, = 1 A^. Then normalized obligors' stock returns have N- 

dimensional normal distribution w.th zero means, standard deviations equal to I and 

correlation matr.x (4 4) Denote by ,J ,heir pint density f^nct.on. Then the 

joint probability that at time /, obligors will be in credit classes R,^Ji,^^ R^^ respectively 

is given by 

^{K='h <=^^...K = /<. < = /^J 

''v;,-.('..'>)\;.^.,(/,./,) 

''i.yif'-'.) V:.„('»./,) Kl.f^Wj.) 

(3.6) 

4. Stochastic processes for index and stock returns under Q measure 

In this section we derive stochastic differential equations for the index and stock returns 
under the martingale probability measure O. 

Let be an instantaneous risk-free interest rate at time /. We assume that is a 
determm.stic positive function of Mhen ,ndex values T^. - = 1,2 A/, and stock prices 

"^''^ ^' for non-dividend-paying stocks (this assumption can be relaxed) satisfy 
following equations; 
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dl =V"'(r,d( + d^.dh] + 0-:,dB^ + . . . -r d"j,dBf ) , / > 0 , 

(4.1) 

f;,"' = r"', w = I,2 V/, 

dH"; = //; (/;<:// + o-,", + fT.",c//^; + . . . + o-,"„ J/^," ) , / > 0 , 

(4.2) 

H;; = h", ii = \.2 ,v. 



B;' = B' + jr\ls, c/ = I /;, (4 3) 

are Brownian motions under the probability measure (J defined by the Radon-Nikodym 
derivative: 

f ==xpf-2j,w-ijKf.Al. 



|r,f.(r;Mr;r*... + (.,"r. 
where vector-function r, =(r'.r," rj'J solves equations: 

^ "! ^' + ^"'2 +■■■ + 07 1} r,' '=/',"' ~ ,"' = 1,2 M , 

(4.5) 



+ + • • • + o", ,,r, = a' - /; , /; = 



(4.6) 



In the matrix form: 
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cr,r, =a, 

(4.7) 

where fi,=(^i]-r, //,"-/-J. S,=(a]~r, a^' ~r,) . We assume that solution r, of 

(4 7) exists for all / e [/„,/,] 

Similar to (1 5), (1.8) one computes index and stock log-returns from (4 1), (4,2) for the time 
interval [l^,t^]=[V,Jf„\. 

X, = c), + jO^JB, .v,_ = 0 , 

(4.8) 

y,-f>, + jo;c/B, ,/„</</,,>', = 0 , 

(4.9) 

where 

«.'k «,"r 

dB^={dB\ dB^) . 

Let Fy[) cov<?{.,.) be an expected value and covariance matrix under the probability measure 
Q, conditional to -algebra . Then from (4.8), (4.9) follows that 
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(4.10) 



i^'iy.hf^. i4.u) 

cov^' (.V, . .V, ) = cox'ix, .x,)= ]efi]ds . (4 12) 

cov<^{y, .x,) = coviy, . -v, ) = i o-^^ V.V , (413) 

cov<-'{y„y,) = cov{y„y, ) = } cr^crfc/v . (4 14) 



In particular. 



var^'(;-;)=var(>-,'') (4 j5) 

for all 77 = 1,2 /V . Moreover, 

cov<^{£,.x,) = cov{6;.x,) = 0 , (4 16) 

cov^'{£,.6;) = cov{£,.£,), (4 17) 

but 

/^''k) = ^,-M-/?>0. (4.18) 



5. Credit migration model under Q measure 

In section 3 we defined the relationship between credit migration events and movements of 
the normalized stock returns. This allows to compute probabilities of credit migrations for 
different credit classes under the O measure. 
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Let obligor .V" e S be in the credit class = at time /„ , According to the assumption in 
section 3 its credit rating at time will be R" = if and only if u-," e {Kj{'nJi\h"^.,{t,,j^)\, 

y = l'2 K , where thresholds //"^(/o,/, ) are derived from the standard normal distribution. 

Then probability under O measure for obligors .S",^' S"' to migrate from the credit 

classes Rl = /< /^^ = at time /„ to the credit classes 7?,'^ = R^^ R'^' = R^ at time /, 

respectively is 

q(rI = R„ < = R,^ \Rl = /e„ R;i = R,^ ) 

(5.1) 

As follows from the definition of if, = (u' w^'J and (4. 15) 

(5.2) 

where 

(5.3) 

and r," (/„,/) is a constant given by 
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Then for the vector = (u,' .j; / using (3 4), (4 14), (4. 15) one has 

cov'^()r, ,»■,)= cov'^(h,,)i J 

= /->'(z;',Cr,-At>',))(>'-/':{v,).yz,:) 
= z,:',cov^'(j„^,)z,.', 
= z;' cov(>',>',)z;', 

= COv(u', ). 

(5.4) 

From equations (4.9), (5.3), (5 4) it follows that the random vector u; = (»•' .l^'f has a 

normal distribution under the probability measure (J with zero vector of means, covariance 

matrix cov^' (..■,...,) = cov(w, ,»■,) and probability density function /,(.r,..r, x,) Then 

w" is a standard normal variable for each // = 1 N and 

^c(A,^o=(.X<e(/c(A,./,)./7:,,(/,,,/,)j) 

=cXv^>r,:(/,,.Oe(/;:,(/,,0./',v,(/,,,ol) 

= Qk e {Kju,)-cr,{'MK, ,(/n,/,)-c:(/,„/,){ ^^^^ 

= M^:.-,(/,>.0- M^:. (/,../,)- r:(/,„/,)). 

Therefore, similar to transition probabilities ^ /,), transition probabilities ^ /,) 

under martingale measure () are defined by shifted thresholds 
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8:.M,.',)= K. (/.../>) 



(5.6) 



Note that transition probabilities ./."^(z., ^ /,) depend on /, and can be different for different 
obligors Equations (5.5), (5.6) show the relationship that exists between two transition 



matrices > 



Similar to (3.6) for the O probability of joint migrations we have 

q(K=^ K = \.K=^ K = 

=c'fe:-fe..('„.o-^;,,,.,(/o.ol....<eU- (5 7^ 



6. Construction of transition matrix under 6* -measure over risk horizon 



In this section we suggest, as an application of the theory developed in earlier sections, a 
practical way to obtain a O transition matrix for the period [/,„/,], t,=t^ + H, where H 
denotes a risk horizon, typically ! year. We suppose that we are provided with a P transition 
matrix and a set of O default probabilities, one for each initial credit rating, over the period. 
The matrix obtained this way will be used in subsequent sections. 
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The main point here is that of implementation issue. In practice, the differences between 
default thresholds of the given P transition matrix and those determined by the given 0 
default probabilities may not be independent of the credit ratings, contrary to the theory (see 
(5.6)) This is the case usually because the P transition matrix is obtained empirically from 
past historical data while the O default probabilities are obtained from current price data. 
Below, we make a modification to the given P transition matrix and then apply the theory to 
obtain a O transition matrix over the risk horizon Note that the last column of such O 
transition matrix will coincide with the given (J default probabilities. 

First, we introduce some general notation. For obligor .V" eS and observation times />.y, 
.V,/ e T, let ({"{s i) denote the (J default probability of .V" over the period [.v,/] given its 
initial rating of at time .v , i.e., 

V," (v ^ /) = Q(It" = D I R: =R,), / = 1,2 A' . 

Recall that the default state /^^ =■ D is absorbing so that <y;;.(.y -> /) = 1 . Put ic = A'-l, and 
package them in a k -dimensional vector 

V(v^/)' 

we have q"(.v -> /) > 0 , where 0 denotes the k -dimensional zero vector. 

Now suppose we are given a P transition matrix T," = ' A"; !' for each individual obligor 
S\ respectively, and O default probabilities q"(/„ /, ) over the risk interval [/„,/,]. Put 
= (/"('n /, ) and r," =N -'{(/;•), 1 < / < k . Let 
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as in Section 5 Define transition matrix = || p'-^ ) || 

jN(z: +(':), j = K 



/V,(/n,/,) = 

J<K (6.1) 
where o, .iHiXl^^ P'U,,J^)-0 if , < A" and /^;:. (/.„0 = 1 . We have adjusted 
the last column of the given .atnx /; ■ so that the new r default thresholds satisfy 
'(//■,.(/,„0)=-.r,: for ,<A-, ie.. they differ from the O default thresholds .-; by 
constant T,'; The other entries corresponding to the probabilities of migration into non- 
default states are then adjusted while maintaining, for each row, the same relative weights 
with respect to the survival probability, i e.. 



^- !>:,:{'. jy \- JK; '-'^^ (6.2) 

Now. we appiv the theory of change of measure to to obtain the desired O transition 



matrix over the risk horizon 



With /•; - II (.,„.,) II of (6 .), let ) be its thresholds. In particular, we have 

A;(/„./,) = ^-'(A",(/„,/.)) = .-;+r,: 

(6.3) 

Let = II (/,„/, ) II be the corresponding transition matrix under Q ; by (5.5) and (5.6), its 
(/,./) -th entries are defined by 

<M,J^) = N{gl ,(/,„/,))-A^(^,"^ (/,„/,)), /, ^,,2 fc ^^^^ 
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where = . Using (6.3), (6.4) and '(-.) = 0 , we have 

' ^ construction, the last column of fits exactly the given 

default probabilities under (J over the period (/,„/,]: 









N/:;. (/,„/,) 












^ 1 







(6.5) 



In sun,.ary, a su„abl. adjus,.=m ,s .ado ,o .he given /- ,ra„s,„on n,a,n.v sueh ,ha, „s la. 
column is compatible w„h ,he given O default probab,l„ies under ,he iheory, , e,. ,he,r 
thresholds di to by a eons.an, C independcn, of i„„ia, c,ed„ raung Cass For eaeh row, ,he 
remaining emnes in ,he /■ ,rans„ion ma,r,x are also adjusted by rerain.ng same relative 
weights with respect ,„ the survival probability The constam C is then subtracted from other 
non-default thresholds of the modiHed /• matrix to obtain the corresponding thresholds, and 
therefore the entries, of the desired (J matrix 

Note, that ,„ most praccal applications transition matrices for individual obl.gors are no, 
available Instead only one common for all obligors /' transition matrix is g,ve„ I, can be, for 
instance, transition matrix published by one of rattng agences l,ke Moody s or S&P, or the 
internal bank's transitton matrix. For such case following modification of the procure 
above can be suggested We assume tba, O default probabilities q-„„ , are also obligor 
independent so that ,"(,„ ^ )= 0 Le, be the , -th component of „,„ and 

put 
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Here, for C'„, we may want to average the differences --, over only the /-th terms for 
which 7, and p, j; are not too small (so that the thresholds are away from - k); also, instead 
of taking the average, one can choose ('„ which gives the least mean square error. Similar to 
(6.1), define the new transition matrix T'' =\\ p, j('„Jt)\\ by 



iN{z,+C„), j = K 



where u. = ' /^.j ('»'',) = 0 if J<K and /\-.;^- (/n.'i ) = i Therefore, the last 

column of the given matrix 1], is adjusted in such way that the new 1' default thresholds 
satisfy A^"' (A ^- (/„,/,)) = ^, + C„, i.e., they differ from the O default thresholds z, by 
constant C„ independent of i . Similar to (6 2) we have 

P..j('»J^) ^ P,.j ^^^^jj iJ<K. 
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7. Forward transition matrices under ( 



In .his s«,„„, we sugges, a„ algori.hm ,ha, builds a senes of ,ra„si,i„„ .a.nces ™,chi„, 



lout this 



.he B,ve„ ,e™ s.^cure of defaul, probabilities over various ma,ur„ies Through, 
«c.io„. we assume ,ha, probabili,ies and .ransi.ion matrices are under O-^easure Such 
series will be useful, for e«.ple, when pricing credi, derivatives via ris. neu.ral valua.ion. 

Recall .ha, for each obligor .V . S we assu.e ,ha. i.s credi, ™gra,ion process observed a, 
.imes in f for.s a MarKov chain on s,a,es R under y -^asure. We view f as 'universaP 
in ,ha, all observa,ion u.es considered fro. now on are ele.en,s of f ,n pracice, we have 
a panicular fi„i,e subse, of observauon „mes known from beginning. 

Le, T = {„„ <„, <...<„J (cf)bease,ofobserva,ion,imes. By a ,,/,™,,„„„ 
™<™<,.v „„ T , we mean a sequence of , ransi.ion ma,rices over ,he periods 

[».,.«,).(«,.«,). .(„„„„„j 
Wi.h such series, .he ,ra„s,.,on ma,rix over any period („„.„,, w„h „<* is gi„e„ by .he 
produc.of,ransi„onma,rices over [,,,.,,,., J K ,.„,). respecing , he order 

Le. T. = = 0,1.2 L ) be a se, of 'sundard' observa,io„ ,imes, where H is ,he 

risk horizon, we assume ,ha, T, = f For obligor ,V.S, suppose we are given a term 
s.mc.ure of defaul. probabili.,es under ^-measure. By appropria,ely ,„,erpola,i„g if 
necessary. ,his means .ha. we are given ,he O defaul, probabiliiies sa.isfying 
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where e is the k -dimensional vector whose components are all equal to I . For h , 0 < b < L 
and an invertible transition matrix over the interval K,/J define vectors q'' -^,,) as 
follows: 

M Potent.ally, these vectors may have negative components because 

A/ ' can have negative entries. However, if these vectors satisfy 

0 = q.;', /.J < q'' (/, /,.,)<... £ q; ) < e, (7 2b) 

then we refer to them as the ,„ (foruarJ) dcfaul, [>rnhahilines oh.ainedfrnm the iniUal 
ten. snncno. nf Jefaul, pmhahUmes (7 i) ..v.,, uansu.on nuun. Af Note that the equality 
in the first relation of (7 2b) is equivalent to the condition that the first k components of the 
last column of equals q"(,„ -^/J, i.e., it matches the initially given default probabilities 
for the period [/,,,/J . 



We wish to construct certain series of (J transition matrices 

on T„ such that the following proposition is satisfied with Z, ^ ,->•/) 



Proposinon: Given initially a term structure of time /„ default probabilities satisfying (7 1 ), 



ime 



>"ng 

there is a set of invertible matrices Z,. Z,, .... Z, such that for each A . 0 < A < /, , the t, 

forward default probabilities obtained from (7.1) using M = Z,Z, -Z, satisfy (7 2b), i.e., 
they are non-negative and increasing. 
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Proof of Proposii ion. We shall produce Z,'s recursively At each stage, we aim to keep the 
entries corresponding to the probabilities of non-default migrations as 'close', in terms of 
maintaining their relative weights (This can be replaced by any other reasonable preference. 
For example, instead of seeking to retain same relative weights amongst the non-default 
probabilities, one can skew them in the direction of default probabilities: if default probability 
has increased then one can increase the relative weights towards lower credit classes, etc.), to 
those of the transition matri.x which was produced from the given data for [/„,/,] in 
Section 6 



the given default probabilities for the interval. If Y, is invertible and condition (7.2b) already 
holds for A/ = then we set /;,«-'(/„ ^ /,) = >; . More generally, consider the following 1- 
parameter family of matrices 

Z,(a) = a:\', +(i-a)};, 0<a<l, (7.4a) 

where 



.'•J q"(/., ^/,) 



(7.4b) 

0 1 

with 7, =r/;(/„ ^/j< I and /;' = l-7,"(/„ ->/,) for / = 1,2 k. The matrices Z,(a) all 

have the same last column and in particular is independent of a and equal to that of \\ 



Choose the "smallest" a , call it or,, in (7.4a) such that vectors q;',(/, ^/,) defined by (7.2a) 
with invertible A/ =Z,(q:) satisfy (7 2b) with h = \ there (in practice, we fix a positive 
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integer/>, from beginning, and restrict values of a in the set a(/)=(«,,/.2 
Then is chosen to be the smallest . i. for which condition (7 2b) holds.) That such 
a, ex.sts is clear s.nce for a = | , the condition is satisfied. Thus, we have the t.me , forward 



default probabilities 

»Sq-(', -*':)Sq-(,, ..<,-(,, ^,_)^^ 

■'('.->',) =Z,(Q,). 



(7.5) 

with respect to M = Z,(a, ); here, subseript M is omitted We set 7y(/„ 



Now. for subsequent penods [,,.,.,), „ , , „ 

.ransitton matnees. Using the t,« forward default probabtltttes (7,5, as given, we ea„ 

obtain in a similar way an invertible transition matnx Z, on [,,.„ whose last column is 

compatible with the (forward, default .probability ,-(r, ^„ „„„ 

■n each row. the same relative wetghts as those of ,or. as before, have the wetgbts 

distributed according to some preference, Note that the time , forward default probabilities 

obtatned from the time r, default probabilities ,7 5, using of,,,,,) w.ll be the same as 

.he ones obtained from the ttme ,,, default probabilities (7 1, using Z.Z, „„ [,,„,^,, 

Repeating the procedure for intervals f/, / 1 \, / i 

I'-' J, I/,, /J and so on, we construct a desired series 

of transition matrices, (7.3), on T„ . 



The key point about the construction is that the series of transition matrices built up this way 
fits the g.ven term stn^cture of default probabilities over each maturity 0 < / . , , in that 
the cumulative trans.t.on matnx Z,Z,..Z, on the interval has (the first . 

components oO the last column matching the given q''(/„ -^/,), 
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Remark There are many choices for such series We have suggested simply one method that 
is practical and reasonable in incorporating much information about the P transition matrix 
given for the initial period as well as preserving non-negativity and increasing of the forward 
default probabilities defined by the constructed transition matrices. Other alternate methods 
for producing such series can be made available. Actually, in practice, one may prefer to be 
given as exogenous input a series of forward O transition matrices in place of the initial term 
structure of () default probabilities 



8. Joint process under (J on T niid sipplicntioii to pricing of credit swaps 

In the previous section, we have built series of Q transition matrices on standard observation 
times T„ for single processes R" . In application, for example in pricing of credit derivatives, 
we need to 'translate" this series to cover a different set of observation times 
T = {'/„<//,<■•■< //„} (T includes nonstandard' time steps corresponding to payment 
dates, times to instrument's maturity, etc., which need not belong to T„) for a inint process 
For this, it is sufficient to describe it for single processes, for then the joint migration Markov 
process will be determined by (5 7) on each period ^'J as the co-migration of the 
obligors are governed by the correlation among their asset returns. We make the assumption, 
in lack of any better information, that this correlation is constant in all periods, but this too 
can be given as input and made different for different periods. For the single process, it 
essentially requires that we provide a way to 'interpolate' a transition matrix for a given 
period to define a transition matrix for a subinterval of that period. We shall simply take a 
certain Taylor series expansion of the 'fractional' power of the transition matrix, the power 
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being proportional to the subinten^al's length, and when necessary make a minor adjustment 



in some cases. 



Let Z be a transition matrix for the period [/,./.,.,] and suppose [«, r] e [/.„/,., ] ,s a proper 
subinterval; put 9 = -,J Now consider the matrix Z" defined by 

where ((6)). = 9(9 - , )..,0 - A- . ,)/., as in the Tayior expansion around the identity matrix 
/ When senes (8 ,) converges, we may tn^ncate it after a certain number of terms to 
evaluate with a small error. Note that a typical transition matrix, especially over a period 
that ,s not too long, i.s 'near' / In actual implementation of (8.1), it may give nse to 
possibly negative entr.cs, though small m s.ze. In th.s case, one makes an adjustment in the 
resulting matnx, for example by replacing such entries by 0 and decreasing the adjusted 
amount from the diagonal probability. In most cases, this procedure works well and is robust. 

We summarize various constructions made thus far Stan.ng with a transition matrix on 
the risk horizon and an initial term structure of O default probabilities we construct series of 
transition matrices for single processes on standard observation times. These are then 
translated to produce series on a prescribed set of obse,.at,on times T, first for the single 
processes and then for the joint process. Of course, for single processes, the O senes thus 
obtained can be used to imply the corresponding P series if differences in the.r thresholds 
are provided. 



Appiicanon: We indicate how the constructed series of O trans.tion matrices can be used in 
evaluating credit swaps. A credit swap on a corporate bor d or a reference facility provides 
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protection to the holder of the reference bond from the loss of value when the bond issuer 
defaults. The buyer of the contract makes a series of premium payments, contingent upon the 
survival of either the seller or the bond issuer. Here, the scheduled dates for the premium 
payments are used as observation times T and the joint process will be that of the asset 
returns of two obligors, namely the counterparty C and the issuer B of the reference facility, 

The price of a credit swap can be computed as the expected value, under fj -measure, of the 
difference between fixed rate and contingent payments. Also, when valuing the credit swap 
for risk measurement, it is important to consider whether the counterparty is a seller or a 
buyer, as the two situations are not exactly the opposite of each other. The payoff function for 
the contingent payments can be given generically in the form 

^rr=Zf.f,: (8.2) 
the sum is over an exhaustive set {rc} of all possible, mutually exclusive events, is the 
indicator function of event /•;, and is the cashflow, usually of what was promised in the 
contract or some fraction of it due to seller's own default, when event /: occurs Similarly, 
one has the payoff function A',,^ for the fixed payments side (buyer) 

Denoting by t,. and t„ , respectively, the random default times of the counterparty C 
(seller) and the reference party fi , then typically E is an event with a specified range for 
values of t^. and t^, . For example, 

E = ( event that u, , < < //, and r^. > ii, ) 
where //,.., < n, are two specified times, usually a pair o.^ consecutive premium payment 
times. Depending on £■, is a function involving notional amount or forward values of the 
reference facility, a strike price, recovery rates of B and :', etc. and depends on whether 
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according to the credit swap's specification. 

The pHce „r.He co„,„,e« Cai. U .he e.pec,e« aisc„„„,e. cas.„„.. 
pricing measure O . It is equal to 

Given a ,e™ ..„„e „r ,„,e.c. ™e. .e „. .ee disc.™,i„,, , ^ ^„ 

app.„.™,e. .eas„„a.v .n a s,..,,,f„™a.d .a„„er. Thcefore, ..e ™a,„i„, .ey 

even-s per:ai„,„, ,„ ,He defauU .i.es ,. and . ,„ p,a„i.e, a„ ,.e p,„,a.i,i,es 

.rans,„o„ ™,r,ces. We .„„s,ra,e ,h,s, for example, for ,he following even,: 

= ( event thai „„ <r, s,,.., and < ) 
^e. /„ = («..«,. ) be the initial ratings of («,0 at time „„ and «, the rating of . at U.e 
"' ^™''>'ly for C. , etc) Then, we have 

ai') = C(S, ,=/J.C, =/J|/„) 

'Zt'r. =A«.., =/J,(-., 

P-e.es Note that the proba.llt.es appearing ,„ the „na, stage are entHes of the transition 
™.nces. over per,ods „„..,.„. , and ,„,„..,, „«ce O,., are easil, a.ta.ned fro. the 
.ransitio„.atHees.Thisco.p,etesthe„na,pi.eofeo.pntatio„,„ pricing credit swaps 
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Part 3 - Stochastic Rates 
1. Introduction 

Credit VaR model with stochastic interest and FX rates (CreditVaR II) is an extension of the 
Credit VaR model described in the Part I (referred fiirther as CreditVaR I). The chart in Figure 
6 describes the implementation of CreditVaR II model. CreditVaR II is a purely simulation 
engine and relies on generating correlated risk factors without pre-generation (like CreditVaR 
I), and fiill re-pricing of the portfolio for each scenario. Outputs are Value-at-Risk and 
various other statistics for portfolio as a whole as well as standalone and marginal statistics 
for each instrument in the portfolio. 
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2.1. Simulation of Base Curves 

Let 0 be the current time and h be the risk horizon. Let >■(/,/') be the continuously 
compounded zero coupon rate seen at time / for a bond maturing at l+T. Note that T is the 
term to maturity, not an absolute time point. Let f{iJ\J\ ) be the forward rate seen at time 
/ for a contract starting at / + 'l\ and maturing at / + /;,. We assume that the zero curve 
maturity buckets are the same for all currencies and all credit ratings, and we denote the 

maturity terms by « = {/,,'; /J, where /„ = 0 < /, </, <;■••</,. For convenience, we 

denote the zero curve as seen at time /; by y = (>',,>• ) = iyihj,),v(h,t, ),-■■, v(h,t )) 

and the one seen at time zero by y" = (>-,", >\ >■;') = (KO, /,),>'(0, /,),•••, v(0,/ J) The 

zero rate for a term between two bucket points is determined by linear interpolation (If / < /, 
or / then the zero rate is set to v, or \\ , respectively.) 

Our first objective is to generate samples of y, which are correlated within itself, spreads, FX 
rates and equity returns. (Note that there is a vector y for each currency We dropped the 
subscript for currency to avoid heavy notation.) Consider the forward rates seen at h, 
denoted byx= (.r,,x, -vj =(/(/»,/,„/,),/(//,/,,/,), -,/(/;,/,.,, /J) . x can be given by 

x, = and jr, = ~-^'-'^ .. ! . for / > 1 . (u 
'.-'.-\ 

Conversely, 

-^.=fl^^(^ (2) 

Written in matrix forms, 
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(3) 



Denoting the above transition niatrix by T, then x' = T 'v' where th • 

y where the inverse matrix is 

given exphcitly as 



(4) 

The foHowin, method is proposed to generate interest rate scenarios. Let x" = . • 

= (m.,.),/(0,,,,,..,,,0,,.,,,.theWd rates seen at time 0^ 

- Wd rates x are ,ognorma„y distributed and that the ,og relative changes 

'^"'"^ a m.t.variate norma, d.stnb.ion 

"^^^ ^ = - is the covariance matrix. There.re, the 

forward rates are guaranteed to be positive. 
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= (,r"exp^ -^ft'exp^J, then aoply equation (I) to get continuously compounded zero 

rates y 

2.2. Simulation of Spreads 

Given a credit rating /e . / = I A' - 1 , let .v, be the spread curve for R, over that of the 

base curve .v Thus. >- + .v is the zero curve for /e rating The restrictions on variability of 
the spreads arc the following: a) the corporate curves should increase with decreasing the 
credit rating; b) the corporate curves should not overlap for any time point / We propose the 
following scheme for generating stochastic credit spreads In order to reduce the number of 
risk factors in the simulation and satisfy the restrictions above, we generate the spread curves 
by a combination of Monte Carlo simulation and linear interpolation The risk factors 
involved in simulation of spreads are: 

- Credit spread .v, (/, ) for the highest credit ratmg for (typically Aaa) and for maturity 
equal to the first bucket point on the zero curve (/, ) 

- incremental spreads Av, (/, ) = .v, (/, ) - .v ^ (/^ ) , , = 2 K - I for maturity /, 

■ ^A. -i('*) = ^\ i('J- -^/. -i(0. ' e„ the difference between the spreads for the lowest 
credit rating for the longest and shortest maturity. 

The simulation relies on the lognormal assumption for the above risk factors After they are 
generated as a part of a simulation scenario, we obtain spread curves by the following 

procedure: 
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^ 0.,ai„ = , ^^„_^, ^ ^^^^^^^^^^^ ^^^^^^ 

are incroising when credilradns is decreasing. i,e, .v,(,,)>.v, ,(,,) fo,/ = 3..,y, 
2 0.,ai„ .,e ^^^^^^^^ ^^^^^^^^ 

for any scenario 

3. Obtain .he res, of ,he spreads for .a,un.y . by partitioning ,he in,e„a, [o.,, ,,,;] 
proponionaiivro ,(,,, by .a.ing ^.O^^'.O.J^ This way we ge, spreads .ha, 
increase when credi. ra.ing declines for ma.nri.y „. 

4. ""'ain. he in.erior Of, he spread curves by linear i„,erpo,a.io„be.ween ,(,,, and .,,(,.,, 

The corpora.e.ero curves a. Hs.hori.o„ are obrained by adding .He approp„a.espread.o,He 
base curves 

2.3. Simulation of FX rates 

The FX rates at risk horizon are assumed to be iognormaiiy distributed. 
2.4. Simulation of credit migrations 



credi, n„gra„ons and defan„s are sinrulared s,n,i,ar,y ,o ,he n,e,hodoiogy of Credi,VaR , 
(see Pa. i). The assun,p,ion ,s ,ha, credi, n,igra,ions are generared by ,he change in asse, 
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value of the obligor, which is linked to industry participation and idiosyncratic component. In 
the framework of CreditVaR II, the migrations are based on the transition matrix, 
idiosyncratic obligor returns and the correlation of industry indices with the rest of the risk 
factors and within themselves. 

2.5. Statistics 

Statistics contain means and standard deviations for some risk factors (forward rates, spreads, 
FX rates) and the covariance matrix of all the risk factors including industry indices. 

3. Input Dntn 

Obligor data: 

Information about obligors is organized into a database containing details of their credit 
ratings, industries, and countries. 

Purlfolio Jala: 

Information about financial positions is organized into portfolios of exposures. It covers 
different types of instruments such as fixed income instruments, loans, commitments, letters 
of credit, etc. 

Market data and transition probabilities: 

These include yield curves, spread curves, foreign exchange rates, transition probabilities 
from one credit rating to another for different credit rating systems (Moody 8 states. Moody 
18 states, S&P 8 states) 
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S.a.i.ics co„,.„ .ea„. an. .a„d„d .ev.a.io. r„. f^ccs (f„™.,. .... spread. 

FX rates) and ,hc covariance .aTix of all ,he risk facors including indusTy indices. 



4. Risk Statistics and Measurement 

Credi.VaRllp.„ducesriskn,easn.n,en,s .Ha, .es™,la„c ..cones produces .yCdnVaR 
1^ On .l,e portfolio level, ,hc following quantities arc calculated: 

- Present value and Forward value of the portfolio at risk horizon 

■ Ois,r,hu.,o„ of the ponfoho value at risk horizon and its .can. standard deviation, 
skewness and kurtosis 

Value-at.Risk nuntbers at various percentile levels with 95% confidence bands 

■ Expected Excession (or Expected Loss Given Default, of the percentile level 

On the instrument level, the following is calculated and reported: 

• Present value and Forward value of the ins.mmen. „ „sk honzon 

- Mean and Standard Deviation of the tnstruntent value at nsk horizon 

■ Marginal Standard Deviation and Delta-Standard Deviation (see below, of the 
instmment for a given portfolio 

- Capital amount (see below) 

- Optionally: Standalone Value-at-Risk numbers at various percentile levels 

- Marginal Standard Deviation of instniment 

- Delta-Standard Deviation of instrument 



SUBSTITUTE SHEET (RULE 26) 



wo 01/52121 



79 



PCT/CAOO/01440 



5. Credit Migration Model 

Credit quality migrations are modelled according to Credit VaR I methodology ( Part 1). 

6. Simulation Engine 

First, the engine prepares the necessary data for simulation. Threshold levels of standardized 
asset returns representing credit rating changes are determined for each obligor in the 
portfolio using the transition probabilities; this part depends only on the obligor s credit class 
and not on the obligor itself Relevant data is read from the database and is stored in the 
appropriate data structures in memory. 

The scenarios are generated based on the covariance matrix and other statistics of the risk 
factors that are relevant for a particular portfolio. Each scenario is composed of the base 
curve part, spread curve part, FX part and credit migration part. 

For each of these scenarios, portfolio has to be revalued. In CreditVaR I model it was done 
quickly using pre-processed prices of instrument. CreditVaR II does full revaluation of the 
portfolio, without cash flow bucketing or any portfolio compression, based on the set of 
generated risk factors. Each time a default occurs in scenario, i.e., when the sampled asset 
return value of an obligor is below the default threshold level, a random recovery rate is 
generated according to a beta-distribution whose defining parameters are governed by 
seniority associated to the obligor's instruments. Finally, we obtain a distribution of portfolio 
values and from it the relevant risk statistics. 
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7. Implementation 

The Global Analytics Credit VaR M program utilizes the methodology described above. The 
core computational engine is implemented in C++ Currently, interface for Windows is being 
developed, using the ADO database technology Windows version currently utilizes a 
Microsoft Access database, ADO technology allows CreditVaR II to be used with the 
majority of available databases 

It is noted that those skilled in the an will appreciate that various modifications of detail may 
be made to the preferred embodiments described throughout this specification, which 
modifications would come within the spirit and scope of the invention as defined in the 
following claims. 
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We claim: 
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' A .e.h„d or an adj„s.e. 

migration comprising the steps of 

• obtammg a transition matrix under probability measure (/'), 

• adjusting tbe transition matrix to be consistent .i,H deraui, pr„babi,i„es under 
maningaie measure ((,) by making a coiumn of default probabilities in .be 
.ransition matrix consistent with default probabilities calculated under (,;, and 

• scaling the other entries in the adjusted „, transition matnx to compensate for the 
adjustment while retaining the relative weights among non-default classes from the 
obtained (/', transition matrix in the adjusted (/>, transition matrix, 

A credit risk model for the analysis of complex credit instt^ment. comprising: 
an obtained transition matHx under probability measure (P). and 
an adjusted (/, transition matrix based on the obtained a, transition matrix and 
con.ste„t with default probabilities under mea.re (0 by making a column of default 
probabilities in the adjusted ^ transition matrix consistent with default probabilities 
calculated under (0, 

w.ere,n the other entries in the adjusted (.) transition ntatrix are scaled to 

compensate for the adjustment While retaining the relative weights among non-default 
Classes from the obtained (/>, transition matrix in the adjusted (/>, transition matHx 
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3 The credit risk model of claim 2, further comprising a Monte Carlo simulation engine to 
generate valuation scenarios based on probability distribution computed using the adjusted 
(P) transition matrix. 

4 The model of claim 2 or 3, wherein the model is implemented in computer software on a 
computer readable medium. 

5 The model of claim 4, wherein the software is for use on compatible hardware 

6 The model of claim 2 or 3, wherein the adjusted (/') transition matrix is used to price credit 
swaps 

7 The method of claim I, wherein the step of adjusting the obtained (P) transition matrix is 
performed in accordance with the following: 

R = {/^, Rfj \ - are all credit classes 

P,: - corresponds to default state 

= p,.j( - adjusted (/') transition matrix for the time intcn.-al [/,„/,] 

P,.,( J - transition probability from credit class /< to credit class R 




/\.y('..<'i) = 0 '^^ J<K and p^^ii„J,) = \ 
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' "^'''"^'^ P^°b^^'i'ity of default for the credit class R, 
(from the obtained {P) transition matrix) 

" probability of default for the credit class 

■ ^ °^ '^'■^'^'t classes (can be the whole set R = 
or some subset 

- denotes the number of elements in G. if G is empty, put 

C, = 0 

- standard normal cumulative density function. 



N{] 



The method of claim I, where, the step of adjusting the (/>) transition matrix is performed 
accordance with the following; 



« -\\p:J'„J.)\\ - adjusted (/^ transition matrix for the t.me interval [,...,,] for 
obligor // 

- obligor-s n transition probability to migrate from credit class 
/< to credit class R 



V^.P.r j<K 

Pk.jOcJ>) = 0 if j<K and plA'n,t,) = 



SUBSTITUTE SHEET (RULE 26) 



PCT/CAOO/01440 



NO 



- obtained (/>) probability of default for the obi.gor . over 
'■"terval [/„.,, J if it is in the credit class /< at time/,, 

- (Q) probabWty of default for the obligor „ over interval 
[/„./,] if it is in the credit class /< at time /„ 

- instantaneous risk-free interest rate at time / 
-obligor „ stocklog-return over interval [/„,,,] 

- obligor /, drift function of time / 

- standard normal cumulative density function 



A method of pricng a credit swap, comprising the steps of 

• — -r.overariskhor.on[,,,.,Ja. 

probabilities over the same period [/,^,/,|, 

'=',..,.,... .0. „e (<;, .ansMon .art, a„C a„ iniUa, ,er. s.™c,.e of 
default probabilities under O over [/„,/,] 
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• Translating the first series to produce a second series, having premium payment dates 
as the observation times by working with 'fractional' powers of transition matrices for 
subintervals; 

• Deriving a third series of transition matrices from the second series for co-migration 
on credit classes for pairs of obligors; 

• Multiplying cash flow for each payment date by a risk-free discount factor and by 
probability contained in the third series and summing the computed series for all 
payment dates; and 

• Subtracting the sum from a payoff fijnction for fixed payments, 

10 A method of pricing a credit swap, comprising the steps of 

• Obtaining a (P) transition matrix over a risk horizon [/,„/,] and the (0 default 
probabilities over the same period [/,„/,]; 

• Creating an adjusted {P) transition matrix and (0 transition matrix over [/,„/,]; 

• Building a first series of forward (0 transition matrices for standard observation 
times / = /,./;,/„... from the (0 transition matrix and an initial term structure of 
default probabilities under O over [/u,'i]; 

• Deriving a third series of transition matrices from the first series for co-migration on 
credit classes for pairs of obligors; 

• Translating the third series to produce a second series, having premium payment dates 
as the observation times by working with 'fractional' powers of transition matrices for 
subintervals; 
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• Multiplying cash flow for each payment date by a risk-free discount factor and by 
probability contained in the second series and summing the computed series for all 
payment dates; and 

• Subtracting the sum from a payoff function for fixed payments, 

1 1 . A method of pricing a credit swap, comprising the steps of 

• Obtaining a (/') transition matrix over a risk horizon and the {()) default 
probabilities over the same period 

• Creating an adjusted (P) transition matrix and (0 transition matrix over [t.,J^], 

• Building a second series of forward (0 transition matrices for premium payment 
dates from the (0 transition matrix and an initial term structure of default 
probabilities under O over [^„/|]; 

• Deriving a third series of transition matrices from the second series for co-migration 
on credit classes for pairs of obligors; 

• Multiplying cash flow for each payment date by a risk-free discount factor and by 
probability contained in the third series and summing the computed series for all 
payment dates; and 

• Subtracting the sum from a payoff function for fixed payments, 

12, A system for pricing a credit swap, comprising: 

• Means for obtaining a (/') transition matrix over a risk horizon and the (0 
default probabilities over the same period [/(,,/,]; 
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. Means fo, creating an a<,«.d .™s,.io„ 

. Means f„, f.o. ,He ad,ns.ed .ansUion n,a,„. and an ,„i,a, ,er. s.™a„. „r 
defanl, prCabUl.ies under „ „ve, (,,,). , „„, 

matrices for standard observation times / = /,./,. 
• Means To. .a„s,a.,„, ,He „„, se.es .„ produce a second series, having pren^iu. 
paymen, da.es as .he observa„on .taes by working .i., 
transition matrices for subintervais; 
. Means for deriving a third ser.es of transition matrices from the second series for co- 
migration on credit classes for pairs of obligors, 
. Means for discounting cash now for each payment date over the period by a risk-free 
discount factor and probability contained in the third series, and summing the 
computed series for all payment dates; and 
. Means for subtract.ng the sum from a payofT function for fixed payments. 

A system for pr.cing a credit swap, comprising the steps of 

Means for obtaining a (/, transition matrix over a risk horizon f,,,J and the (0 
default probabilities over the same period [r„j,]- 

Means for creating an adjusted (/, transition matrix and (0 transition matrix over 

Means for bu.lding a Hrst series of forward (0 trans.t.on matr.ces for standard 
Observation times ^^^^ 

structure of default probabilities under O over [/„,/,]; 
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• Means for deriving a third series of transition matrices from the first series for co- 
migration on credit classes for pairs of obligors; 

• Means for translating the third series to produce a second series, having premium 
payment dates as the observation times by working with 'fractional' powers of 
transition matrices for subintervals; 

• Means for multiplying cash flow for each payment date by a risk-free discount factor 
and by probability contained in the second series and summing the computed series 
for ail payment dates; and 

• Means for subtracting the sum from a payoff function for fixed payments. 



A system for pricing a credit swap, comprising the steps of 

• Means for obtaining a (/') transition matrix over a risk horizon [/,„/,] and the (0 
default probabilities over the same period [/„,/,], 

• Means for creating an adjusted (P) transition matrix and (O) transition matrix over 

• Means for building a second scries of forward {(J) transition matrices for premium 
payment dates from the {(J) transition matrix and an initial term structure of default 
probabilities under O over [/,„/,]; 

• Means for deriving a third series of transition matrices from the second series for co- 
migration on credit classes for pairs of obligors; 

• Means for multiplying cash flow for each payment date by a risk-free discount factor 
and by probability contained in the third series and summing the computed series for 
all payment dates; and 
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• Means fo, subtracng ,he sun, fron, a payoff function for fixed pay„,en,s 

instructions to. 

• Ob.a,n a ^ ^^^^^^ ^ ^ ^^^^ ^^^^^^^ 
probabilities over the same period [/„./,], 

. Create an adjusted (/>) transition matnx and ,0 transition ntatrix over [,„.,,]. 

• Buiid a first sertes of fonvard (O, tr^sifion tnatrtces for standard o.servauon ti.es 

"'" """" - .ern, structure of default 

probabilities under (J over [/„,/,] 

• Translate the first senes to produce a second senes. hav.ng pretniun, pay„e„, dates as 
* obsen-ation tt.es by workin, wtth .fractional- powers of transit.on matrices for 

subintervals, 

• Dertve a third ser.es of ,rans,t.on matrices fro. the second ser.es for co-n,i,ration on 
credit classes for pairs of obligors; 

probabihty contatned ,n the thtrd series and su. the computed senes for al, payment 
dates; and 

• Subtract the sum from a payoff fcnction for fixed payments, 

.6 A method of modellm. stochast.c .nterest rates, credtt spreads and fore.n .change rates 
comprising the steps of: 

• Computing obligors- credit ratings at ris. hon.n .ased on counts/industry market 
'ndex logreturns and idiosyncratic components; 
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• Computing base zero interest rates for all currencies and given maturities /„,/ / 

from base forward rates /('„.0./(V:)---/('/,-i.'J corresponding to time intervals 

k./,]. [v,.^]; 

• Computing credit spreads .v,(/,) = .v,. ,(/,) + A.v,(/,) for maturity and credit classes 
/ - 2....K - I (by currency), 

• Computing credit spread ) = s^,^(l^) + ^s^_^(l^ ) for the lowest credit rating 
and longest maturity (by currency); 

• Computing the rest of the spreads for maturity by panitioning the interval 

[O.-V \('>,)] proponionally to .v,(/,) by taking -v/^; ^ .v//, /^^" (by currency); 

• Computing spreads for other maturities by linear or any other type of 

interpolation between •s-,(/,) and s^{l^) (by currency); 

• Computing the corporate zero curves at risk horizon by adding the appropriate spread 
to the base curves (by currency); and 

• Computing foreign exchange rates. 

17. A credit risk model for the analysis of portfolios of credit instruments, comprising: 

• Means for computing obligors' credit ratings at risk horizon based on 
country/industry market index logretums and idiosyncratic components; 

• Means for computing base zero interest rates for all currencies and given maturities 

^M.A from base forward rates /('o.O-/('i.^:) A't-i^'h) corresponding to time 

intervals ['„,A|], [a,,/,] k.i.'J; 
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• Means for computing croUt spreads .v (/,) = 5 r/ ) + aw/ w 

A',; + for maturity /, and 

credit classes/- 2... X-l ;by currency); 

• Means for computrng credit spread ...J.J = s^J, j . , 
credit ratmg and longest maturity (by currency); 

. Means for comput.ng the rest of the spreads for matur.ty , by pan.t.on.ng the 



.Oj = s.Oj^^ (by 



interval [O,s,,0j] proportionally to .v (,,) by tak.ng . 
currency); 

• Means for computing spreads for other matur.ties by l.near or any other 

type ofinterpolation between .•(,,) and s,(,J (by currency). 

• Means for computmg the corporate zero cu.es at ris. horizon by adding the 
appropriate spread to the base curves (by currency); and 

• Means for computing foreign exchange rates. 



-es. ™*e, i„d. ...„s, idiosyncratic co.po„=„,3, ,..e a ,„,„, p™,,w„, ,3„,„„„„ ^ 
.h=. cedi, migration process and in.eres, and foreign exchange ra.es are correlared, 

'9. THe credi, r.s. .ode, of Car. ,S, a,nher con,pr,s,„, a Mon.e Cario ..„,a„„n engine ,o 
generate valuation scenarios based on Monre Carlo s,n,u,a.,„ns of risk factors. 
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Obligor Data 

• Credit ratings 

• l\/larket participation 

• Idiosyncratic data 



iVIarket Data 
Yield curves 
Transition probabilities 
Exchange rates, etc. 




Portfolio preprocessing 
Pre-computing exposure values 
Compute asset correlations between the 
obligors 

Aggregation by obligor and by seniority 



Analytic Engine 

• Compute mean and standard deviation of 
exposure forward values tor a given time 

• Compute marginal standard deviation and 
marginal risk by obligors 

• Compute total standard deviation and total 
mean of portfolio forward value for a given 
time honzon 



Simulation Engine 
• Compute mean and standard deviation of 
exposure forward values for a given time 

honzon 

Compute total standard deviation and total 
mean ot portfolio fora/ard value for a given 
time horizon 

Compute Value-at-Risk at various percentile 
levels together with confidence bands 




Final Report 

• Report the results 

• Mean and standard deviation obtained by 
two different methods should be close 



Figure 1 
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Obligor Data 

• Credit ratings 

• Market participation 
' Idiosyncratic data 



Market Data 

• Spot base curves 

' Spot spreads 

' Spot FX rates 
Transition probabilities 
Market data statistics 
(mean + covariance matrix 
of base curves, spreads, FX 
rates and industry indices) 



Portfolio Data 
• Bonds, Loans 
Floating Rate Notes 
Interest Rate Swaps 
Loan commitments 




Simulation Engine 

* Sr/fl!- '^°"''"'*"9 °' correlated credit 

rnr^f 9;^"°"= (^"d defaults), corporate zero 

Re-pnce portfolio for each generated scenario 
■ distribution of portfolio value at time 

?^S".'^!h"- Value-at-Risk 
InS T^'^^r'^ ''^"•^s). Expected Excession of 
portfolio value at time honzon for g.ven percenSle 

Ss't^len^ ^"^ "^^^9*"^' ^^'^'^^ -ch 



Final Report 
• Customized reports of the simulation results 



Figure 6 
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